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A  requirement  to  determine  analytically  the  acceleration- 
time  relationship*  for  the  cargo  compartments  of  CV-2 
and  CV-7  aircraft,  when  subjected  to  various  crash 
conditions,  and  to  determine  the  effect  of  these  rela¬ 
tionships  on  a  load-limited  restraint  system  provides 
the  basis  for  this  report.  The  techniques  -used  In  this 
report  represent  a  logical  extension  of  previous,  similar 
programs;  In  the  absence  of  test-crash  data,  the  results 
contained  herein  may  be  used  with  a  high  degree  of 
confidence. 

The  excessive  load- limiter  displacements  (2  to  16  feet) 
Indicated  for  certain  crash  conditions  considered  In 
this  report  are  Judged  to  be  operationally  unacceptable 
by  this  activity.  However,  the  load- limiter  concept 
represents  a  significant  safety  and  operational  improve¬ 
ment  over  the  conventional  restraint  devices  for  most  of 
the  crash  conditions  considered.  This  activity  suggests 
that  a  trade-off  study  be  conducted  to  determine  the 
maximum  safety  that  can  be  achieved  by  load- limited 
restraints  commensurate  with  acceptable  operational 
factors. 

Since  the  Army  no  longer  has  operational  responsibility 
for  the  CV-2  and  the  CV-7  aircraft,  this  activity  has 
not  planned  further  effort  in  the  area  of  this  report. 
However,  results  of  this  contract  will  be  forwarded  to 
the  appropriate  Air  Force  agency  for  consideration. 
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SUMMARY 
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This  report  presents  the  findings  of  an  investigation  into  the  crash 
pulse  of  fixed-wing  cargo  aircraft  and  the  resulting  behavior  of 
cargo  restrained  by  load  limiters. 


A  crash  pulse  simulator  computer  program  was  developed  that  ob¬ 
tains  acceleration-time  histories  at  selected  stations  in  the  cargo 
compartment  and  under  various  crash  conditions.  This  simulator 
was  employed  to  obtain  crash  pulses  for  a  wide  range  of  input 
parameters,  both  for  the  CV-2  and  the  CV-7  Army  aircraft.  The 
resulting  acceleration  pulses  were  studied  to  determine  a  suitable 
spectrum  of  realistic  pulses. 


The  crash  pulse  program  was  subsequently  modified  to  include  a 
routine  that  would  simulate  cargo  dynamic  behavior  during  the 
crash  sequence,  employing  the  floor  acceleration  data  as  it  is 
developed.  This  latter  program  was  applied  to  CV-2  and  CV-7 
aircraft,  under  significant  crash  conditions,  to  obtain  the  dynamic 
response  of  cargo  to  the  crash  pulse,  f 
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INTRODUCTION 


This  report  represents  a  continutation  of  the  analytical  investigation  of 
acceleration-4me  pulses  for  the  CV-2  and  CV-7  aircraft  as  described 
in  USAAVLABS  Technical  Report  65-30,  "Cargo  Restraint  Systems  for 
Crash  Resistance".  The  crash  impact  parameters  and  limits  thereof 
used  in  this  program  were  felt  to  represent  the  conditions  most  likely 
to  be  encountered  in  an  accident  for  the  aircraft  involved.  In  addition, 
it  was  felt  that  all  of  these  conditions  represented  a  survivable  accident. 
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CRASH  PULSE  SIMULATOR 


DESCRIPTION 


The  computer  program  to  simulate  dynamic  behavior  of  a  fixed-wing  air¬ 
craft  accepts  as  input  (1)  accident  configuration  variables,  (2)  structural 
data  for  the  specific  aircraft,  and  (3)  soil  behavior  parameters.  The  ac¬ 
cident  configuration  variables  consist  of  velocity,  sink  rate,  impact  angle, 
and  angular  velocity  (if  present).  The  structural  input  data  include  flex¬ 
ural  and  axial  stiffness  properties,  mass  distribution  data,  buckling, 
yield  and  failure  criteria,  and  various  geometric  data.  The  soil  reaction 
parameters  relate  the  interaction  forces  between  aircraft  and  ground  to 
the  depth  of  penetration,  the  rate  of  penetration,  and  the  horizontal 
velocity. 


The  simulator  output  is  in  the  form  of  computer-plotted  acceleration¬ 
time  curves  for  three  stations  in  the  cargo  compartment.  Additionally, 
the  simulator  displays  the  maximum  plowing  depth,  the  groove  length, 
the  maximum  values  of  generalized  coordinates  for  vibration  modes  ex¬ 
cited  by  crash  forces,  the  maximum  stresses  in  the  fuselage,  and  the 
kinematic  data  at  the  completion  of  rebound  (or  slide-out). 


ANALYTICAL  BASIS 


The  simulator  operates  in  either  of  two  general  modes:  the  elastic  mode 
or  the  plastic-hinge  mode. 


The  elastic  mode  is  assumed  if  fuselage  (  longitudinal)  stresses  are  below 
buckling  or  plastic  limits.  In  this  mode  the  dynamic  behavior  of  the  air¬ 
craft  is  divided  into  two  parts:  First,  the  rigid  body  behavior  is  obtained 
from  the  soil- structure  interaction  forces  and  rigid  body  equations  of 
motion.  Second,  the  vibration  modes  excited  both  by  the  soil  reaction 
forces  and  by  the  rigid  body  motion  inertial  forces  are  evaluated. 


The  vibration  amplitudes  are  found  by  means  of  a  normal  mode  analysis 
in  which  the  first  and  second  transverse  and  the  first  and  second  longi¬ 
tudinal  modes  of  fuselage  vibration  are  considered.  The  amplitudes  of 
these  four  vibrational  modes  are  taken  as  generalized  coordinates  in 
Lagrange's  equations  of  motion  for  the  aircraft.  For  the  assumption  that 
the  damping  coefficient  is  proportional  to  mass  distribution  along  the 
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fuselage  (an  assumption  which  introduces  little  error),  the  Lagrange's 
equations  reduce  to  a  set  of  independent  one-degree-of-freedom  equations: 


where 
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A  solution  satisfying  initial  conditions  is  found  in  the  convolution  integral: 

/l 

Q  (T  )exp  I"  -i(t-T  )1  sinu  (t-T  )dT 
n  I  c  I  n 

n  n  w  J 


which  is  evaluated  numerically  in  the  simulator  program  (reference  Ap¬ 
pendix  I).  With  the  generalized  coordinates  evaluated,  the  vibrational  kin¬ 
ematics  are  determined  and  may  be  superimposed  on  the  rigid  body  kin¬ 
ematics  to  obtain  resultant  accelerations.  The  maximum  bending  stresses 
are  determined  from  the  bending  moments  associated  with  the  transverse 
vibrational  modes  of  the  fuselage. 


The  plastic-hinge  mode  of  aircraft  behavior  is  followed  by  the  simulator, 
if  at  some  fuselage  station  the  effective  buckling  or  the  yield  stress  has 
been  exceeded  (when  operating  in  the  elastic  mode).  In  the  plastic-hinge 
mode,  the  aircraft  is  treated  as  two  rigid  bodies  joined  by  a  plastic  "yield 
hinge"  which  may  transmit  hinge  reactions  as  well  as  a  plastic-hinge  bend¬ 
ing  moment.  In  this  mode,  elastic  vibrations  are  ignored,  as  they  are  con- 
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sidered  to  be  negligible  when  compared  with  the  large  plastic  deformations. 
Again,  the  equations  of  motion  are  solved  by  numerical  integration. 

The  simulation  is  terminated  by  any  of  the  following  conditions: 

1.  The  rebound  has  been  completed;  that  is,  the  soil  penetration 
has  returned  to  zero. 

2.  The  horizontal  component  of  velocity  has  been  reduced  to  zero. 

3.  The  deflection  angle  at  the  plastic  hinge  has  reached  a  critical 
value,  implying  a  fuselage  break. 

4.  The  aircraft  rotates  outside  a  specified  angular  limit  (imply¬ 
ing  an  overturning). 

A  flow  chart  and  an  outline  of  computational  operations  of  the  simulator 
program  are  provided  in  Appendix  I, 
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DEVELOPMENT  OF  INPUT  DATA 


The  input  data  for  crash  simulation  of  the  CV-2  and  CV-7  Army  aircraft 
are  divided  into  three  categories:  accident  configuration  variables, 
structural  data,  and  soil  reaction  parameters.  An  example  of  input  data 
for  a  simulation  run  as  developed  for  the  CV-7  aircraft  appears  in  Ap¬ 
pendix  IV. 


ACCIDENT  CONFIGURATION  DATA 


The  accident  configuration  variables  used  were  those  specified  in  the 
contract  statement  of  work. 


STRUCTURAL  DATA 


The  mass  distribution  data  have  been  obtained  from  weight  analyses  made 
by  the  manufacturer.  These  data  were  adjusted  in  accordance  with  var¬ 
ious  cargo  and  fuel  level  conditions  for  several  simulator  runs. 

The  structural  stiffness  data  may  be  conveniently  separated  into  two 
groups:  that  which  deals  with  the  vibrational  behavior  of  the  aircraft, 
and  that  which  relates  imposed  forces  to  structural  deformation. 

Consider,  first,  the  vibrational  behavior.  The  mode  shapes  and  natural 
frequencies  (for  both  longitudinal  and  transverse  vibrations)  have  been 
obtained  by  means  of  a  separate  computer  program  employing  standard 
techniques  (see  Appendix  VI,  Program  VIBRAT).  Flexural  rigidities, 
longitudinal  stiffness,  and  mass  distribution  are  the  input  to  this  auxil¬ 
iary  program;  these  input  data  have  been  obtained  from  analyses  per¬ 
formed  by  the  manufacturer.  The  output  of  the  auxiliary  program  be¬ 
comes  input  for  the  simulator  program  and  consists  of  normal  mode 
shapes,  natural  frequencies,  generalized  masses ,  and,  for  the  trans¬ 
verse  vibration,  normalized  bending  stresses. 

Next,  consider  the  load  deformation  input  parameters  for  the  structure. 
A  double  modulus  relationship  has  been  postulated  between  crushing  force 
and  structural  deformation  in  the  local  area  of  contact  between  the  ground 
and  the  aircraft.  The  crushing  force  moduli  have  been  obtained  by  com¬ 
puting  average  local  plastic  buckling  loads  for  those  structural  members 
that  fail  by  local  instability  and  by  adding  an  appropriate  percentage  for 
the  resistance  offered  by  bending  of  longitudinal  members.  Based  upon 
observed  behavior  of  a  plastically  deformed  structure,  an  estimate  has 
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been  made  of  plausible  "apringback"  deformation  associated  with  un¬ 
loading.  A  range  of  10-  to  20-  percent  springback  was  postulated  for  a  2- 
to  3-  foot  depth  of  structural  crushing  in  the  region  of  the  underside  of  the 
nose.  The  springback  is  obtained  in  the  simulator  program  by  employing 
a  greater  modulus  for  unloading  as  well  as  for  initial  loading  below  a 
critical  force  (marking  the  onset  of  plastic  crushing). 


Figure  2.  Plot  of  Normal  Force  Versus  Structural  Deformation 

The  onset  of  plastic  crushing  is  signaled  by  a  critical  value  of  applied 
normal  force  (see  Figure  2).  This  critical  value  is  estimated  considering 
those  members  that  would  first  undergo  general  buckling;  however,  this 
value  must  also  be  consistent  with  the  elastic  modulus  for  deformation, 
the  plastic  modulus,  and  the  assumed  points  through  which  the  curves 
must  pass. 

The  simulator  program  computes  compressive  bending  stresses  along  the 
top  of  the  fuselage  for  each  time  increment.  When  a  critical  effective 
buckling  stress  is  reached  at  a  given  station,  a  "yield  hinge"  is  considered 
to  be  formed  in  the  fuselage  and  the  mode  of  simulation  changes.  The 
critical  fuselage  buckling  stress  has  been  computed  on  the  basis  of  the 
critical  stress  for  a  longitudinally  stiffened  curved  panel.  The  computed 
value  has  been  increased  by  an  estimated  factor  to  account  for  the  dynamic 
overload  capability  associated  with  short-duration  pulses  and  redistri¬ 
bution  of  stress  as  a  plastic  hinge  is  formed. 

Once  a  yield  hinge  forms,  a  plastic -hinge  moment  offers  substantially 
constant  resistance  to  further  bending  at  the  hinge.  This  plastic-hinge 
moment  has  been  computed  on  the  basis  of  resistance  offered  to  continu¬ 
ation  of  buckling  and  crushing  of  skin  and  longitudinal  members. 
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In  addition  to  the  structural  data  described  above,  geometric  input  defining 
the  longitudinal  contour  of  the  underside  of  the  fuselage  is  required.  These 
data  were  obtained  from  scale  drawings  provided  by  the  manufacturer. 


GROUND  PARAMETERS 


It  should  be  noted  that  an  accident  can  occur  in  a  wide  variety  of  soils 
whose  properties  may  vary  considerably.  Consequently,  the  ground 
parameters  were  varied  over  practical  ranges  for  a  number  of  simulator 
runs. 


The  normal  soil  reaction  force,  Fjj,  was  considered  to  be  composed  of 
several  components.  First,  the  essentially  elastic  behavior  provides  a 
contribution  that  increases  with  penetration.  However,  as  both  surface 
contact  area  and  ground  resistance  per  unit  area  are  roughly  proportional 
to  depth  of  penetration,  the  assumed  relationship  is: 


where 


N 1 


■  C1ZG 


Z^.  =  ground  penetration 
Cj  =  soil  elastic  modulus 


The  constant  C  ^ 


has  been  estimated  for  a  typical  soil. 


A  second  contribution  to  normal  force  is  associated  with  the  phenomenon 
of  planing  (hydroplaning)  is  associated  both  with  horizontal  velocity  and 
depth  of  soil  penetration.  The  relationship  is: 


N2 


=  C.,Z  X 
2  G 


where 


X  =  horizontal  component  of  aircraft 
velocity 

=  planing  parameter 

The  planing  coefficient  Cg  has  been  roughly  estimated  for  various  soils, 
based  upon  momentum  exchange  and  soil  compressive  resistance.  The 
degree  to  which  the  nose  of  the  aircraft  remains  a  "plane"  as  opposed  to 
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a  "scoop"  would  further  affect  this  parameter.  In  the  computer  simu¬ 
lation  runs,  a  range  of  practical  values  was  used  for  this  parameter. 


The  third  contribution  to  the  normal  force  component  is  that  of  the  direct 
momentum  exchange  effect.  As  the  aircraft  penetrates  the  ground,  it 
accelerates  soil  mass  to  its  velocity.  A  momentum  exchange  relationship 
leads  to 


N3 


‘  2  2 
■  C3ZGZG 


where 


Z_  *  time  rate  of  soil  penetration 
G 

C.j  =  soil  impact  parameter 

The  parameter  is  approximately  equal  to  the  mass  density  of  the  soil 
multiplied  by  the  contact  area  and  divided  by  the  penetration  Z  , 

The  total  normal  component  is  the  sum  of  the  three  separate  contri  - 
but  ions: 


+ 


The  tangential  component  of  the  soil  reaction  force  is  similarly  com¬ 
posed  of  separate  contributions.  The  first  of  these  is  simple  friction. 


where 


f  =  coefficient  of  friction 


To  the  extent  that  the  soil  has  been  penetrated  and  some  manner  of  scoop 
has  been  formed  at  the  nose  of  the  aircraft,  two  additional  contributions 
to  the  tangential  force  would  exist.  The  first  of  these  is  associated  with 
soil  "drag"  or  "plowing"  action  and  is  proportional  to  the  horizontal 
velocity  and  ground  penetration: 


C4XZG 


The  plowing  coefficient  for  a  given  soil  can  only  be  determined  ex¬ 
perimentally,  For  this  purpose  a  simple  experiment  was  designed  and 
conducted  to  provide  a  rough  separation  of  the  plowing  and  friction  phe¬ 
nomena  with  two  extremes  of  scoop  conditions.  The  results  of  this 
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experiment  were  employed  to  provide  order-of-magnitude  values  for  the 
soil  reaction  parameters  f,  C^,  and  C^.  The  experiment  is  described  in 
Appendix  III,  The  other  contribution  associated  directly  with  the  scoop 
effect  is  the  horizontal  momentum  exchange.  It  would  thus  vary  with  soil 
penetration  and  with  the  velocity  squared;  that  is, 

FT3  =  C5X2ZG 

The  coefficient  C  would  depend  upon  effective  scoop  area  and  soil  den- 
sity  and  has  been  roughly  computed.  In  simulation  runs,  this  parameter 
will  be  varied  considerably  to  simulate  conditions  ranging  from  a  de¬ 
finite  scoop  in  a  freshly  plowed  field  to  the  other  extreme  of  planing  over 
hardpan  soil  or  a  concrete  ramp.  The  total  tangential  force  is  then  the 
sum  of  the  individual  contributions: 

FT  =  FT1  +  FT2  +  FT3 


i 
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APPLICATION  OF  CP  ASH  PULSE  SIMULATOR 


The  CV-2  and  CV-7  aircraft  crash  pulses  were  developed  by  the  simu¬ 
lator  program  for  combination^  of  the  following  parameters:  velocity, 
sink  rate,  attitude  angle  at  impact,  soil  conditions,  and  aircraft  weight 
conditions. 

1.  Velocity  was  varied  from  SO  to  120  feet  per  second. 

2.  Sink  rate  was  varied  from  10  to  30  feet  per  second. 

3.  Attitude  angle  was  varied  from  3  to  15  degrees. 

4.  Two  basic  soil  conditions  were  considered:  hardpan,  and  a 
soft  soil  equivalent  to  a  cultivated  field. 

5.  For  each  aircraft,  two  weight  conditions  were  considered: 
operational  light  and  operational  heavy. 

A  summary  of  significant  results  is  contained  in  the  following  tables, 
which  relate  maximum  cargo  compartment  acceleration  (excluding  short- 
duration  peaks),  time  duration  of  pulse,  and  velocity  change  to  the  var¬ 
ious  input  parameters. 

An  example  input  to  the  simulator  program  and  developed  output  is  con¬ 
tained  in  Appendix  IV. 


TABLE  I 

RESPONSE  TO  IMPACT  VELOCITY 


Velocity 
(ft/ sec) 

Soil  Condition 

Maximum 

Acceleration 

(G) 

Pulse 

Duration 

(sec) 

Velocity 

Change 

(ft/sec) 

80 

Hardpan 

6.  1 

.  138 

12 

100 

Hardpan 

6.  4 

.  140 

12 

120 

Hardpan 

7.  0 

.  144 

13 

80 

Soft  Soil 

9.  4 

.  188 

26 

100 

Soft  Soil 

11.4 

.  200 

33 

120 

Soft  Soil 

12.  7 

.  220 

42 

CV-7  Aircraft,  operational  light. 

Sink  rate  =  20  fps;  impact  angle  =  12  degrees 
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TABLE  II 

RESPONSE  TO  SINK  RATE 


Maximum 

Pulae 

Velocity 

Sink  Rate 

Soil  Condition 

Acceleration 

Duration 

Change 

(ft/sec) 

(G) 

(aec) 

(ft/ aec) 

10 

Hardpan 

3.5 

.  126 

7 

20 

Hardpan 

6.4 

.  140 

12 

30 

Hardpan 

9.  1 

.266 

26 

10 

Soft  Soil 

6.  9 

.  248 

26 

20 

Soft  Soil 

li.4 

.  200 

33 

30 

Soft  Soil 

14.4 

.268 

54 

CV-7  Aircraft, 

operational  light. 

Velocity  *  100  fps;  impact  angle 

=  12  degrees 

TABLE  III 


RESPONSE 

TO  IMPACT  ANGLE 

Impact 

Maximum 

Pulse 

Velocity 

Angle 

Soil  Condition 

Acceleration 

Duration 

Change 

(deg) 

(G) 

(sec) 

(ft/sec) 

3 

Hardpan 

6.  1 

.  126 

11 

9 

Hardpan 

6.2 

.  124 

11 

12 

Hardpan 

6.4 

12 

15 

Hardpan 

7.2 

14 

3 

Soft  Soil 

9.0 

.  144 

19 

9 

Soft  Soil 

10.  6 

26 

12 

Soft  Soil 

11.4 

.  200 

33 

15 

Soft  Soil 

12.  2 

.244 

41  . 

CV-7  Aircraft,  operational  light 
Velocity  *  100  fps;  sink  rate  = 

20  fps 

l 
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CARGO  SIMULATOR 


DESCRIPTION 


The  computer  program  to  simulate  cargo  restraint  systems  is  designed 
to  obtain  the  dynamic  response  of  a  cargo  retention  system  (employing) 
load  limiters)  to  an  applied  acceleration  pulse.  The  cargo  simulation  is 
accomplished  not  as  a  separate  computer  program  but  rather  as  a  sub¬ 
routine  appended  to  the  crash  pulse  simulator  program.  As  the 
specific  crash  pulse  is  generated,  it  becomes  input  (internally)  to  the 
cargo  simulation  subroutine. 


MATHEMATICAL  MODEL 


The  cargo  is  assumed  to  rest  on  the  cabin  floor  (with  a  suitable  friction 
coefficient  between  cargo  and  floor),  restrained  fore  and  aft  by  load 
limiters  as  illustrated  in  Figure  3. 


Figure  3.  Cargo  Restraint  by  Load  Limiters 

Each  load  limiter  provides  unidirectional  restraint;  the  aft  load  limiter, 
for  example,  offers  restraint  to  forward  displacement  of  the  cargo  rela¬ 
tive  to  the  floor  (when  slack  ie  taken  up). 
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The  longitudinal  force  applied  to  the  cargo  is  transmitted  from  the  air¬ 
craft  through  either  the  aft  or  the  forward  load  limiter,  as  well  as  through 
the  cargo  floor  by  means  of  friction.  The  maximum  value  that  this  force 
may  assume  is  the  load-limiter  limit  force  plus  the  coefficient  of  friction 
times  the  normal  force  between  floor  and  cargo.  (The  normal  floor  force, 
in  turn,  depends  upon  the  normal  component  of  floor  acceleration.) 

At  each  time  increment,  the  subroutine  accepts  as  input  the  generated  ac¬ 
celeration  components  (longitudinal  and  normal)  of  the  floor  at  the  cargo 
location.  For  the  case  of  zero  relative  velocity  (of  cargo  to  floor),  the 
cargo  mass  times  the  longitudinal  acceleration  is  compared  with  the  com¬ 
puted  maximum  possible  longitudinal  force.  If  the  former  is  less  than  the 
computed  limiting  value,  no  relative  acceleration  occurs.  Otherwise, 
the  cargo  acceleration  is  the  computed  limiting  longitudinal  force  divided 
by  cargo  mass.  The  resulting  relative  acceleration  is  then  the  difference 
between  the  floor  acceleration  and  the  absolute  cargo  acceleration. 

For  the  case  in  which  the  cargo  has  a  relative  velocity,  the  applied  cargo 
force  becomes  the  load-limiter  load  plus  the  computed  friction  force 
(assuming  tie-down  slack  is  not  present).  Again,  the  applied  cargo  force 
determines  the  absolute  cargo  acceleration  and  in  turn  the  relative  ac¬ 
celeration. 

Numerical  integration  of  the  kinematic  relationships  is  used  to  obtain 
relative  cargo  velocity  and  displacement  (at  each  time  increment)  from 
the  above -computed  relative  acceleration.  The  load-limiter  stroke  is 
finally  obtained  as  the  maximum  cargo  displacement. 

A  flow  chart  and  computational  relationships  for  the  cargo  simulation  sub¬ 
routine  are  given  in  Appendix  II. 
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APPLICATION  OF  CARGO  SIMULATION  PROGRAM 


The  cargo  simulation  program  was  employed  for  ten  representative 
crash  pulses  (as  shown  in  Table  IV).  For  each  pulse,  the  load>limiter 
stroke  was  computed  for  five  different  load-limiter  limit  forces,  cor¬ 
responding  to  2G,  4G,  6G,  8G,  and  10G  limiting  accelerations.  The  re¬ 
sulting  limiter  strokes  are  tabulated  in  Table  IV  along  with  the  input  con 
ditions  for  each  simulation  run.  (The  center  of  the  cargo  compartmeht 
was  selected  as  a  representative  location. ) 


A  family  of  curves  showing  load-limiting  force  versus  cargo  displace¬ 
ment  appears  in  Figure  4. 


Appendix  V  shows  typical  output  sheets  from  the  cargo  simulation  pro¬ 
gram. 
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TABLE  IV 

RESULTS  OF  CARGO  SIMULATION 


Aircraft 
Configur¬ 
ation 
and  Soil 
Condition 

Velocity 

(ft/aec) 

Sink 

Rate 

(ft/aec) 

Impact 

Angle 

(deg) 

Max. 

C 

(C) 

Pulae 

Duration 

(aec) 

Velocity 

Change 

(ft/aec) 

Load 

Limiter 

(C) 

Cargo 

Dlapl. 

(ft) 

120 

30 

15 

20 

.214 

60 

2 

16.  67 

4 

7.  66 

*WD 

6 

4.22 

J 

8 

2.  41 

_ 

10 

1.40 

4 

c  ~ 

90 

20 

9 

11.  6 

.  158 

26 

2 

2.  56 

.2  o 
~  t r 

4 

0.  91 

2  c 

6 

0.29 

8  in 

8 

0.  08 

0 

10 

0.  00 

i 

80 

10 

6 

5.6? 

.  128 

12 

2 

0.  43 

> 

4 

0.05 

6 

0.00 

8 

0.  00 

10 

0.  00 

100 

30 

15 

7.4 

.  300 

25 

2 

1.33 

4 

0.  33 

6 

0.05 

4 

8 

0.00 

•I 

10 

0.00 

120 

30 

12 

7.  5 

.  300 

27 

2 

1.28 

0  1 

4 

0.  32 

*5  -p 

6 

0.03 

u  4 

8 

0.00 

8 1 

10 

0.00 

0 

120 

30 

6 

6.0 

.240 

21 

2 

0.68 

*7 

4 

0.  13 

> 

u 

6 

0.00 

8 

0.00 

10 

0.00 

120 

30 

12 

14.2 

.  300 

64 

2 

13.  10 

4 

5.01 

* 

6 

2.  04 

4 

8 

0.  78 

i 

10 

0.21 

4  £ 

120 

30 

9 

13,2 

.  198 

37 

2 

5.20 

s  0 
.2  « 

4 

2.02 

*■»  *-» 

6 

0.84 

*«  0 

V  V) 

8 

0.29 

£ 

10 

0.08 

r» 

120 

30 

6 

l’.O 

.242 

44 

2 

4.67 

i 

> 

4 

1.66 

o 

6 

0.55 

8 

0.  12 

10 

0.01 

p* 

4 

c 

80 

30 

3 

8.  6 

.  112 

20 

2 

0.  58 

5  e 

4 

0.18 

ts  *  i 

M  >  T> 

6 

0.02 

'Ml. 

>  6-  *  * 

8 

uoxx 

— 

===== 

10 

15 
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Legend 

Weight 

Condition 

Soil 

Condition 

Velocity 

(ft/sec) 

Sink 

Rate 
(ft/ sec) 

Impact 

Angle 

(deg) 

(1) 

Heavy 

Soft  Soil 

120 

30 

12 

<2) 

Heavy 

Soft  Soil 

120 

30 

9 

(3) 

Heavy 

Soft  Soil 

120 

30 

6 

(4) 

Heavy 

Hardpan 

100 

30 

15 

(5) 

Heavy 

Hardpan 

120 

30 

6 

(6) 

Light 

Soft  Soil 

120 

30 

15 

(7) 

Light 

Soft  Soil 

90 

20 

9 

(8) 

Heavy 

Hardpan 

80 

30 

3 

Figure  4.  Plots  of  Load-Limiter  Force  Versus  Cargo  Displacement 
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RESULTS 


The  investigation  has  led  to  the  development  of  a  crash  pulse  computer 
simulator  (Appendix  IV)  providing  an  analytical  tool  for  the  study  of  crash 
pulse  time  histories  of  fixed-wing  aircraft.  This  tool  was  applied  to  the 
CV-7  aircraft  to  obtain  tabulated  acceleration  data  (Tables  I,  II,  III)  for 
variations  of  pertinent  input  parameters,  including  soil  conditions  and 
accident  configurations. 

In  addition,  the  investigation  has  led  to  the  development  of  a  cargo  simu¬ 
lator  program  (Appendix  V)  to  serve  as  an  analytical  tool  to  study  cargo 
acceleration  levels  and  displacements  during  a  crash  sequence.  This 
simulator  was  applied  to  obtain  tabular  data  and  plots  (Table  IV  and 
Figure  4)  relating  cargo  displacements  to  load-limiter  values  for  various 
crash  conditions. 

It  should  be  pointed  out  that  all  of  the  acceleration-time  pulses  generated 
in  this  program  fall  within  the  human  tolerance  limits. 


APPENDIX  I 


CRASH  PULSE  SIMULATOR 
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COMPUTATIONAL  OPERATIONS  (Refer  to  List  of  Symbols) 


1.  Geometric,  structural,  and  soil  parameters  that  depend  upon  the 
angle  0  are  computed  as  polynomials: 

2  3 

H  =  A  +  A  ft  +  A  0  +  A  0 

n  nl  n2  n3  n4 


2.  The  coordinates  x  and  y  are  found: 


1 


x  =.  H2  sin  0  -  (sB  -  —  )  cos  0 

y  =  (sB  -  -pj-  )  sin  6  +  H2  cos  ft 
1  • 

3.  The  interference,  Z,  and  its  time  derivative,  Z,  are  found; 

Z  =  Y  -  y  ,  Z  =  AZ/At 


4.  The  normal  force,  is  computed: 

a.  For  positive  Z  and  F^T  less  than  critical  F^. 


N 

FN  =  K1  <Z  -  ZG> 


where 


B2  -B 


g  =  ci  +  <zG)  c2 


G 

B 


&  ZG/At 

"-SC)H8*  *  K1 


and 


F.,  critical 
N 


1  +  80 
H, 


K, 


1  +  40 
H, 


C.  = 


1+40 
H  . 
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C2  H. 


S  =  scoop  factor  (varying  from  0  to  1) 


b.  For  F_,  >  F_,  critical,  Z  positive, 
N  N 


AFn  =  kKj  (  AZ  -  A ZQ) 


where 


AZ. 


AZ 


[2gZQ  +  (1  -  Sc)  HgX  ]  kKj  4  1 


k  =  reduction  factor  for  modulus  K 


1 


c.  For  Z  negative, 


afn  =  KjUZ-  azg) 

where 

zr 

azg  =  A Z-2 


5.  The  tangential  force  is  computed: 

ft  =  fFN  ♦  sc[h8x  +  h9>c2]zg 

where 


f  =  effective  coefficient  of  friction 
6.  Accelerations  are  computed: 

*  =  FN/mA  •  l  -  FT/mA 
»  -  1  Fnx  t  FT(y*-^))  /lA 

where 

m  .  =  aircraft  mass 

A 

I  .  =  aircraft  mass  moment  of  inertia 

A 
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7.  Generalized  forces  for  vibration  modes  are  computed: 


d  4 


1 


P*1,SJ»  +  M— Vi  W 
s  1 


<*♦ 


1 


P*Z(8J)  +  M_dT(SJ>'  ?  aL  Wi  ♦z(Si) 


where 


c*3(sj>i  ViW 
C*4<8J»i  Vi  W 


P  =  F^  cos  0  -Fj,  sin  fl 

C  =  F  cos  0  +  F  sin  6 
T  N 


M  =  C 


[<v 


+  — )  sin  6  +  x  cos  6 


] 


^  1‘  *  '  ^4  =  norma^zed  vibration  modes  as  functions  of  s 

Sj  =  distance  along  fuselage  axis  to  center  of  force 
application  (from  nose) 

a.,  a.  =  rigid  body  acceleration  components  (transverse  and 
longitudinal  at  station  s. 

u  ,  =  mass  associated  with  i^  finite  subdivision  of  fuse¬ 
lage 


8.  Generalized  coordinates  are  evaluated: 


l  -  r  Qn(t)  *  1  \ 

q(t+At)  =  exp  (-A  At)  jq  (t)  cos  «  it  +  h(t)  +  - . 5—  sin“nAt> 

n  c  (  n  n  L  n  m_  w_  j  f 


n  n 


where 


n  =  varies  from  1  to  4 


S  =  damping  coefficient 

/*L  1 

m  =  generalized  mass  =  /  V4 


dx 
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-  r 

n  ‘  V“n  4 


=  natural  frequency 


and  h  (t>  is  found  from: 
n 


l  r  Qn(t)  A  t 

h  (t  +  a  t)  =  expf-^At)  <  h  (t)  +  —  cos  w  At  -  q  (t)  sinS  At 

n  2lLn  m  w  J  n  n  n 

1  n  n 


9.  The  bending  moment  is  computed: 


d2  ♦. 


mb(8.|  =  0(8.)  [q,— +  q2— T“l)] 

d  8  d  s 


where 


EI(s.)  =  Flexural  rigidity  at  s  =  s. 

10.  The  longitudinal  acceleration  is  computed 

a_.  =  a.  +  a,r. 

Ti  i  Vi 

where 

a^,.  =  total  longitudinal  acceleration  at  s  =  s. 

a.  =  rigid  body  acceleration  at  s  =  s. 

a  =  relative  acceleration  from  longitudinal  vibration 

Vl  2  2 

d  «3  d  *4 

avi  =  q3— 7 r(si>  +  q4— r<’i> 
dt  dt 

11.  Components  of  velocity  and  displacements  are  computed  from  kine¬ 
matic  relationships: 

AX  =  (X  -  ~  )  At 

where 

•• 

,  X  =  updated  acceleration  (that  is,  acceleration  at  end  oi 
time  interval,  At) 


AY,  AO  a?e  found  similarly. 


?.2 


where  X  =  updated  velocity 
AY,  A0  are  found  similarly. 

12.  For  the  "yield-hinge"  mode  of  aircraft  deformation,  rigid  body  accel¬ 
erations  are  computed  as  follows:  The  angular  acceleration,  6p(>>of 
the  forward  portion  of  the  aircraft  and  the  angular  acceleration,  0^, 
of  the  rear  portion  (aft  of  the  yield  hinge)  are  each  computed: 

.  -  G3-  g5.g2.  g6 

F  Den 

Gr  °6~  °3*  °4 
F  ~  Den 


where 


t 


GD  =  DXR(DYF’  eF  4  DYR  ’  6  R  )  ‘  °YR  (DXF’  6  F  +  °YF‘  eR  ) 


I„,I„  =  mass  moments  of  inertia 
F’  R 

m.,m_,m  =  masses 

A  r  R 

X  ,  Y  =  coordinates  of  yield  hinge 
H  H 

M„  =  yield-hinge  mbment 
H 

DXF  -  XF  -  XH 

X  ,  Y  =  coordinates  of  forward  center  of  gravity 
r  F 

DYF’DXR,DYR  =  similar  to  Dxf 

2  2  2 
rF  =  °XF  +  °YF 

2  2 
r„  =  similar  to  r ^ 

R  F 

•  m  • 

13.  Acceleration  components,  X^,,  Yp,  X^,  Y^,  of  forward  and  aft  section 
center  of  gravity  are  computed: 

FT  4  mRC7  mR 


3 

> 

1 

1 

3 

> 

o 

vO 

fn 

•mRC8 

mR 

mA 

mAC10 

ft 

mFG7 

mF 

mA 

'  mAG9 

F 

n 

-mFC8 

mF 

mA 

mA°10 

where 


G7  =  ®F2dXF  +  9R2°XR 
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G8  =  6 F  °YF  +  0 R  DYR 


G9  =  DYF'  0F  +  DYR‘  0  R 


G10  "  °YF  0F  +  °YR  ’  0  R 


14.  Velocity  and  displacement  components  are  updated  from  kinematic 
relationships: 

4Xn 


XF  =  <XF-—  )At 


where 


X  =  updated  acceleration 
•  •  •  •  • 
ayf,  Ae  f,  axr,  ayr.  Ae  R  are  found  similarly. 


aXF  =  (XF-  —  >4t 


where 


X^  =  updated  velocity 

AY_,  A0  AX  ,  A Y  ,  A6q  are  found  similarly. 
F  F  R  JR  K 


! 


APPENDIX  II 
CARGO  SIMULATOR 
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COMPUTATIONAL  OPERATIONS  FOR  CARGO  SUBROUTINE 


1.  Longitudinal  acceleration  is  computed  by  crash  pulse  program. 

2.  Normal  acceleration  (augmented  by  an  equivalent  acceleration  for 
weight)  is  computed: 

••  N  M 

a^„,  =  X  sin  a  +  Y  cos  6  -  0s  ,  +  32.  2  cos® 

FN  J 


where 


s 


J 


Distance  from  center  of  gravity  to  the  cargo  location 


3.  If  cargo  relative  velocity  (relative  to  the  floor)  is  zero  at  beginning 
of  time  interval  A  t,  cargo  acceleration  is  computed  as  follows: 


where 


aF  for  (’FlL  '  fCaFN)<aF<  (RLL  +  tCaFN) 


-FLL  ‘  fCaFN 

for  aF  < 

("FLL  ‘  fCaFN, 

RLL  +  £CaFN 

for  aF  > 

(RLL  +  fCaFN 

a 


C 


a 


F 


R 


LL 


F 


LL 


cargo  acceleration 

floor  acceleration  (longitudinal) 

Floor  acceleration  (normal) 

coefficient  of  friction  between  floor  and  cargo 

rear  load-limiter  limit  load  (in  the  form  of  equiv¬ 
alent  cargo  acceleration) 

forward  load-limiter  limit  load  (equivalent 
acceleration) 


4.  If  the  cargo  has  a  relative  velocity,  cargo  acceleration  is  computed 
as  follows: 


=  R „  .  +  f~a™,  for  v__.  >0  and  no  slack 

C  LL  C  FN  REL 
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aC  =  -FLL  •  fCaFN  tor  VREL-  °  and  "°  Slack 

5.  If  slack  is  present  in  system,  the  cargo  acceleration  is: 


a  =  a 
C  F 


for  -f_a  <  a  _<  f_a„, 
C  FN  F  C  FN 


ac  "  "fCaFN 


for  >p<  .  fc»FN 


a  =  f„a  for  a  >  +  f„a 

C  C  FN  F  C  FN 


6.  The  relative  cargo  acceleration  is  computed  as: 


aREL  “  aF  ‘  aC 


with  a  ^  taken  positive  in  the  for¬ 
ward  direction  but  a^,  and  a^.  positive 
in  the  aft  direction. 


7.  Relative  velocity  and  relative  displacement  are  computed  by  numerical 
integration: 


4vrel  1  (a 


SsREL  =  |V 


A  aREL 

)At 

REL 

2 

avrel 

)  At 

REL 

2 

where 


s  =  relative  (forward)  displacement  of  cargo 

KlL 


v  =  relative  (forward)  velocity  of  cargo 

R  EL 
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APPENDIX  III 


SOIL  DRAG  EXPERIMENT 


An  experiment  designed  to  separate  phenomena  involved  in  the  longitudinal 
resistance  offered  by  the  soil  (during  a  crash)  and  to  establish  relative 
magnitudes  of  significant  parameters  was  conducted  using  a  typical  clay- 
sand  soil. 


DESCRIPTION 


An  instrumented  "shoe"  was  attached  to  the  underside  of  a  heavy  cart. 

The  cart  was  towed  at  various  speeds  over  a  trough  of  prepared  soil,  with 
the  shoe  set  at  a  level  for  a  given  penetration  of  the  soil.  (See  Figure  5.  ) 


Figure  5.  Soil  Drag  Experiment. 


Instrumentation  provided  a  history  of  normal  and  tangential  forces  for 
each  run.  Six  runs  were  accomplished  (and  duplicated);  they  consisted 
of  three  velocities  with  the  shoe  orientated  in  each  direction,  i.  e.  ,  with 
the  sloped  edge  forward  and  the  blunt  edge  forward. 
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The  results  are  tabulated  as  follows: 


TABLE  V 

EXPERIMENTAL  DATA 


•locity 

(fps) 

Shoe 

Orientation 

Tangential 

Force 

(lb) 

Normal 

Force 

(lb) 

4 

Slope  Forward 

550 

800 

22 

Slope  Forward 

520 

700 

38 

Slope  Forward 

750 

750 

4 

Blunt  Forward 

500 

550 

22 

Blunt  Forward 

400 

400 

38 

Blunt  F orward 

450 

350 

If  a  mathematical  model  for  the  tangential  force  is  postulated  as 


=  fF 


N 


v  ASX' 


ranges  for  the  coefficients  f  and  v  may  be  computed  consistent  with  the 
above  data.  These  ranges  are  found  to  be: 


f  varies  from  .  3  to  .  9 

3 

v  varies  from  10  to  46  (lb-sec/ft  ) 
where 

2 

A  is  the  total  contact  surface  (ft  ) 
s 
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INPUT  NOTATION  FOR  COMPUTER  PROGRAMS 


INPUT 

NOTATION 

A( I ,  J)  -  (1=1,  3;  Jsl,  4)  Coefficients  describing  contour  coordinates 
(ft) 

A(I,  J)  -  (1  =  3,  9;  J=l,  4)  Coefficients  describing  ground-aircraft 
interaction 

2 

AS(J)  -  Area  of  structural  cross  section  (ft  ) 

BETA  _  Longitudinal  vibration  damping  factor 

CDTH  -  Critical  angle  of  plastic  hinge  rotation  (rad) 

CIA  -  Mass  moment  of  inertia  of  aircraft  about  "Y"  axis  (slug  ft  ) 

CIM(J)  -  Mass  moment  of  inertia  per  section  (slug  ft^) 

CLN  -  Station  for  cargo  attachment 

CMA  -  Total  mass  of  aircraft  (ft-lb) 

2 

CMU(J)  -  Mass  per  section  (lb  sec  /ft) 

CM(I)  -  (1=1,  4)  Generalized  masses  (lb  sec^/ft) 

CNS.NS  -  Number  of  fuselage  sections 

DF(J)  -  Fuselage  depth  from  longitudinal  axis  (ft) 

DS  -  Length  of  a  section  (ft) 

DT  -  Time  increment  (sec) 


ECA(J) 

ECB(J) 


Normalized  stresses  for  first  and  second  bending  vibration 
modes  (lb/ ft^) 


EC(J)  -  Structural  parameter  (lb/ft) 

EI(J)  -  Flexural  rigidity  of  fuselage  (lb  ft^) 
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EPS 

FLL 

FMU 

MSW 

OM(  J) 

PH1(J) 
PH2(  J) 
PH3(  J) 
PH4(  J) 

RLL 

RPM 

SB 

SNK 

STRE 

THT 

THTD 

VEL 

ZM(  J) 


Error  factor 

Forward  load-limiter  limit  acceleration  (G) 

Friction  coefficient 

Program  switch  code 

Frequencies  of  vibrations  (rad/sec) 

Normalized  deflections  of  first  and  second  vibration  modes 

Rear  load-limiter  limit  acceleration  (G) 

Reduction  factor  to  plastic -hinge  moment 
Location  of  center  of  gravity  (ft) 

Sink  velocity  of  aircraft  (ft/sec) 

Buckling  stress  (lb/ ft^) 

Angle  of  impact  (rad) 

Angular  velocity  (rad/sec) 

Velocity  at  aircraft  (ft/sec) 

Section  moduli  (ft^) 
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APPENDIX  IV 

CRASH  PULSE  SIMULATOR  FORTRAN  LISTING 


PROGRAM  CRASH 

0 1 HFNSION  VClO)<PV(lO)>A(o.4)»H(Q)  .WOO  ( 3Q ) >Q( 20 )  t _ 

1 PQC  70 )  .Phi  (3?)*PH?(*2)*PH3<32)  »PH4(32  I  *CMU(  32  I  »  EC  A  (12  )  »FCR(32) 

- MS  I  Z3.1  »C  I M  <32  )  iAC113+*UL3.±1+B.U.11-U - 

1  OM  (  /■  )  *  C  MUTTs^D  U  )  .CSOC  4l  iCMOM  (4  )  iPLOT(  71  )  .  R  T  (  4  )  .QMX  (  4  1  .  Z  M  (  32  )  ♦ 

1 0  FORMAT ( 8F1 0.6) _ 

RFan  1 0  »  VFL.SNK.ThT.ThTO 

PPAP  1C.  (  ( A (  I  tJI > J=1 >4 1 . T  =  1  .Q  ) _ 

R^AO  1  0  ♦  CM  A«C  I  A.SP  *STRo.COTHt  Rpv.rNS.PS  *FMU»OT 

Rran  10.CMU.CIM.nF.PHl.PH7tPH3.PH4.FCA.FCR.ZM _ 

RFAn  lC.OMtCMfRPTA 

RFAO  300.  ( PLOT (J) »J=1.71 ) .BLNK .PRO  .  AST . EQ.PLS _ 

3jU  FORMAT  C71A1.5A11 

pc/>n  Q.MNTH.NOAV.NYR.NRUN _ 

o  FORMAT  ( 4  I S  ) 

5  C  =  ,  7  S _ 

11  OO  1  2  *  J  = 1  *10 

12  V(J)=0. _ 

V<R)=ThT 

V(4  1=SQRTF(VFL«VFL-SNK»SNK  ) _ 

V(S)=-SNK 

V  (  A  ) -ThTQ _ 

PO  14  1=1.2 

u  V  (  t  i  =  a(T  .1  )+ThT»(A(I  .7)4-ThT»(a(  T.31+THT«a(  1.4)  )  1 _ 

FM-P. 

IF  IThT)  70.71.71 _ 

70  VC  1 1  =  1 ./ (1 ./All.  1  1-600. *THT ) 

71  t i =  sr- 1 . /V ( 1  ) _ 

T?=V( 7 )*S!NF (THT)-T1*C0SFCTHT ) 

VC  2)  =  V(  7  ) *C0SF ( THT  )  +Tl*SINF(  THT  ) _ 

VC  1  >  =  T  7  . . 

FT=oN _ 

MXT  =  1 
NCR=1 
MRTN=1 

THTP  =  THT _ 

OO  17  I =  3 . R 

it  H(n  =  A  C  T  .1  )+THT»(  AC  I  .7  )+THT»  (  A  C  T  .3  )+THT»A  (  1  .4  )  )  ) _ 

H^P  =  H(  3  1  /(  1  .+4,*THT*THT  )  ’  . .  ” 

HAPr ( ] ,+4.  »THT»THT  )  /H(4  ) _ 

HRP  =  M  (  3 ) 

HGP  =  H  C  6 ) / ( 1 ,+0.*THT*THT ) 

H7P  =  H«  8  ) 

H8P  =  H(R  ) 

H1P  =  V( 1  ) 

no  1 3  J  =  1 »?0 

13  0 ( J ) =6 , 

T=n. _ 
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r -  ~ 


wrdmw'-'llij  f< 


- reTrfl-g - 

0>‘C»)  =  OMU)  -  RFTA  *  RFTA/4,0 ___ _ 

0M{4)  =  0M(4)  -  RFTA  *  ofta/4.'-' 

RO  1  F  I  =  1,4 
rtMXU>=0. 

CM O^C I ) = SQRTF ( OM (  I  )  )*CM( I ) 

SND  (  I  )-*ilNF(iO^TF(oM(  I))*nT) 

if  csrx  i )  = cos fm  sqrtfcqmc  n  )»m  > _ 

SWCH=  -1,0 

ro  ->oi  j  =  i  ,ns _ _ _ _ 

»S1  AC(J1  =  0, 

no  i  is  i  =  ,4 _ 

SNO  (II  -  FXPF  (  — R  FT  A  /  ? ,  *  PT  )  *  SNPC!) 

l±  rSQ  (  1  )  =  FXPF(-°FT  A/7,  «  PT)  »  CSQ(I) _ 

PT(1>=0. 

PTC?) *0. 

»Tn  j  srfta 

0T(4 ) ePCJA 

- Trmr. - 

FNr  =  0.0 

fc=«.o 

FU  =  2.0 _ _ _ _ _ 

- 

ZPG  =  0. 

- - 

FTRM  r  0. 

PRINT  ?O0,MNTHtNnAY»NYR,NPUN 
?r0  FORMAT  (  1  Hi  ,  ?0X  «  fiH  PATF  ,7?4,?0X.fh  RUN  *  T  F  /  /  ) 

PRINT  ?0l 

?0i  FORMAT ( fX»48H  TIMF  STA.  A  STA.  R  STA.  C  CACCFL.  G  UNITS)/) 

- MSW  =  1 - 

STRL=0. 

00  11?  1*1.3 

n?  vc  i  )=vn+t  i*nT*v(  i ) 

1 10  PFN=FN 

oo  111  1=1,10 

111  ovc  T )=V( I ) 

THT  =  V ( 3 ) 

DO  1 13  1  =  1  .9 

in  H(  M=A<  I  »1  ) +THT* ( A (  I  ,?)  +TMT#  (  A  (  I  »3  )  +THT»A (  1,4))) 

IF  CTHT)  117.1  18*118 
1  1  7  HC.l  1  =SH 

US  T1=5B-1.Vh(1.) - 

H  (  4  )  =  1  ,  /H  (  4  ) _ • _ - _ 

T?=m(?)«61NF(THT)  -T  1*F0SF(  THT  ) 

H(?)=H(?)*rOSF(THT  I +T 1#S INF  (  THT ) 

H( 1 ) =  T  2 

VC  10)*H<2)-VI?) 

- Zp.WiOi-frVfifl)  >/M - 

Z  =V  c 1 o )  ^ 


1  2  7 

TF(T-.I)  110*700.700 

1  BO 

i  f ( v  u ii  n  a » 1 1 4 » i s  l 

181 

I  F ( Z I  114*114*18? 

1  B? 

fF(vm-.R)  115.700.700 

1  1  4 

PRINT  116 

1  1  6 

FORMAT ( 1 BH  RFflOUNO  COMPLFTF  ) 

on  TO  700 

1  1  « 

SWrw  b  -i.o  *  SWCH 

IF  (bWCH)  118*321*122 

121 

STOP 

l?2 

GA=AC(R)/12.2 

N  =  OA  +  ?l  ,  6 

TF  (Nl  102.102.i0i 

1"  -» 

IF  ( N-7 l 1  104*104*10? 

1^4 

PLOT  (Ml  =  AST 

1-2 

0  R  =  A  C ( 111/12. 2 

<  =  GB  +  26.5 
IF  U)  105.105. i06 


106 

IF  (K-71)  107*107,106 

107 

PLOT  ( K )  =  EO 

10  6 

GO=AC( 181/12.2 

L  =  r,r+i  l  .6 

IF  (L)  10B.108.10Q 

.  109 

IF  ( L-7 1 )  110.110,108 

1  1  0 

PLOT  (L)  =  PLS 

10  8 

PLOT  (26)  «  PRD 

PLOT ( ? 1 JsPRO 

PLOT ( 1 1 1 =PRD 

111 

PRINT  311 *T .GA.OR.GC .PLOT 

FORMAT  (6X,F5.1.3FB.2»71 A1 1 

IF  (N)  112,112.111 

1 1 1 

IF  ( N-7 1 )  114.114,112 

114 

PLOT  (N)  =  BLNK 

112 

IF  ( K )  116,315.116 

116 

IF  (K-71)  117.317.315 

117 

PLOT  (K)  =  BLNK 

1  1  6 

IF  (L)  318,118.110 

1 1  9 

IF  ( L-7 1 )  120.120.118 

120 

318 

PLOT  (L)  =  BLNK 

CONTINUF 

10 

OZ=ZD*OT 

COF  =1,+4.*THT*ThT 

H  <  3 1 SH( 3 ) /COF 

H ( 4 1 =H( 4 ) *00F 

H(6)=H(6)/(2,*C0F-1  .  » 

IFIZG  -  ZGMI  4 1  1 *41 1 *412 

4]? 

ZGM=Zg 

ZSP=Z-ZG 

4U 

pIg  =  Zg 

IF  (DZ)  119.120.120 
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12" 

GO  TO  ( 800,801 ) «  V8W 

fUO 

FUH=H<4  l+Zr>G  »  ZOG/H<  5  1 

r<  =  ( ( 1 ,0-SC )*H(8)*V<4)+1 .0/H( 1 )) /  (  2.0 
ZG*SORTF(Z/  (H(  3  ) »FUH ) +CK*CK ) “CK 

#  FUH) 

zs  =  z  -  zr, 

FN  =  ZS/H( 1  ) 

IF  (FN-1./HI6))  801,803,007 

00? 

MSW  =  7 

801 

fuh  =  h ( 4 1  ♦  Znr,  «  Znr,  /hid 

IF  (FN  -  1./HIM1  804,805,808 

804 

FCH  *  H»3) 

GO  TO  806 

- BTTT 

FCH  *  FC  *  H  (  3  )  “ 

80*, 

OZGsDZ/  (  (  7  .#FUH*ZG+(  l,-Sft*Hlm»VI41l  *FCH+ 1,0  ) 

j5f5o" 

zr,*ZG+nzG 
vs  =  nz  -  dzg 

nFiy  *  dzs/fch 

FN*PFN+DFN 

0010303 

1  1<J 

DZr,  =  0Z*ZG/  (Z-OZ  ) 

FCH  ■  H ( 3  » 

GO  TO  099 

801 

FN8r r ( i,-Sf )*H<0 »«VU)*7.G 

IF  ( FN-FNSC )  808,171,171 

cc 

e 

oc 

FNrFNSC 

121 

FT  3  FHU  *  FN  +  ( H ( 8 )  *  V  C  4  >  +  H(O)  * 

V  C  4 )  *  V ( 4 ) )  *  ZG  »  SC 

Znr,  =  (Zg  -  PZG)/0T 

CTh=COSF(V<3»  | 

STH*S I NF ( V ( 3  )  } 

P*  FN*CTH-FT*STH 

<*  =  FT  *  CTH  +  FN  *  STH 

121 

GO  TO  1 126,1651 ,MRTN 

1  28 

V(71 *-FT/CMA 

PNCR  =  NCR 

V ( 8 ) *  FN/CHA 

V(  o  1  =  (  FN*H  ( 1  )+FT*(V  <2>i-ZG/2.  >  I/CIA 

00  130  1*1,3 

V(i+I  1*  V(  3+1  )  +  (  V  (  f  ♦  I  )  +PV  (  6+  I  )»*OT/2. 

170 

v(n«vi  t  + 1 v ( 3+t  i  +pv i i+n  »*nT/7. 

SJsOS/7. 

DO  131  J*1 »NS 

WOM  J)=V<7 )*STH+V(S )«CTH-V(9)*  (SB-SJ) 

nr 

5J*5 J+05 

DO  132  I  =  1,20 

132 

POI!l«f)III 

$P*SB-( V  ( 7 1 +ZG/2. ) »STH+H ( 1 ) *CTH 

<V<  2 1+ZG/2. >*CTH+H<1  )*STH) 

J«SP/DS+1.0 

om  «  p  *  phi (  j >  +  cmh  *  (phi  (j+i  )-phi  (jp/ns 

0(2)  *  P*  PH2U)  +  CMM  *  (PH2(  J+l  )-PH2(  J)  )/DS 
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0<7)=C*PH7( J) 

0(4  )=C»PH4( J  ) _ 

AOCl  =  V ( 8 (  *  STH  -  V  (  7  )  *  rTH 

SJ  =  SR  -  OS/2* _ 

THSO=V (  6  )  *  V  (  6  ) 

no  137  J  =  I *NS _ 

0(1)  =  Q(1)-CMU(  J)*wnr>(  J|*PH1IJ| 
0(?)  r  Q(?)-CMU(  J)»Wnn( J)»PH?( J) 


0(7)  =  Q(7)  - 

CWU(J)  * 

(  AOFL 

+  THSO  * 

SJ) 

*  PHS(J) 

0(41  =  0(4)  - 

(-MU(J)  * 

(  AOFL 

+  THSO  * 

_SJJ_ 

«  PH4 ( J ) 

i 77  sj  =  sj  -  ns 

00  212  I  =  1.' 

k 

IFIT-DT/2.)  710*210*211 

710  0(4  +  I  )  =  Q  (  1  )  »  r>T  »  OT/(ft«  »  CM  (  I  )  ) _ 

(1(U  +  !)  *  0(4  +  II  *W(5QRTF(0M(  m  *  r>T  ) 

GO  TO  211 _ 

? U  OT  =  PO ( 1 6  +  „H  +  PQ( I ) /OMOMI I)*DT 

0(4  +  1)  =  P0(4  4-  II  *  CSn(I)  ♦  AT  «  SND(I) _ 

o(i a  +  n  *  ot  *  rsod)  -  PQU  +  n  *  SNn(n 
HI  0(  17+1  )  =Q(  1  )/CM(I  )-(0M(  I  )  +RT  (  1  )»RT I  I  )  /  4  « ) *0 ( 4  +  I  l-RT(  1)»PQ(B»1 ) 
T  F  (  ABSF  (0(  4+1  )  )-QMX  i'll)  712*212*214 

214  OMX  (  1  )  =  ARSF  ( Q  (4->  !  _ 

2 1 2  0(fl  +  I)  =  PO  ( 6  (0(17  +  I)  +  P0(12  +  II)  *  DT/2.0 

PSTR=0« _ 

PPSTR=0. 

rMNT  *  CMA  -  TMU( 1 )/7.0 _ 

SJ  =  SR  -  OS/2. 

PSHR  =  0. _ 

pmm=-cmm 

NPP= 1  _ 


401  IF  ( SP-SRfSJ-OS/2. )  407  *407* 40 ? 


407  NPP=2 _ 

SHR=SHR-P 

40?  qMM=BMM-(PSHR+SHR)»OS/?,»CIM( J)»V(9) 
PSHR=SHR 

_ STRsEC4( JI»Q(5)+FCP( J)»Q(ft) _ 

STRN=(7,*STR+6.*PSTR-PPSTR  )  / 8 . 
IF(STRL-STRN1  191*191*190 _ 

1 90  strl=strn 

JSTL»J _ 

191  IF(STRN-STRM)  413*414*414 
414  STRM  *  STRN 

JST  =  J 

413  IF ( STRN  -  STRB )  140*141*141 _ 


rfllMiBMWUWs 


rrr  PFTnrrrrnr~  - 

PSTR^TR 
r-n  TO  144 

141  JM=I-1 _ 

*H  =  ,  40  -  SJ  -04/?.  0 

°RlM  lt?iJH _ 

14?  FORMAT  (?1H  YIFLH  H|Nf,r  AT  STA.  .11) 

MX  T  =  7 _ 

144  rMNT  =  CMNT  -(CMIMJ1  +  rMU(J  +  ]||/’,0 

4^0  5J  =  SJ  -  ns _ 

14*  FORMAT  <1X  *1  on  1  .4/1 

T  »T+OT _ 

r,0  TO  (  1 10.14R)  .MXT 

149  fMF=0. 

T*T-oT 

DO  150  Ns  1  t JH _  _ 

150  C«F*CMF+CMU(N I 
FMp,CMA-rMF 
^LM=C.MR  /CM  A 

TM=rMF/CMA _ 

RMerMR/CMF 

SJ=nS/2. 

SRF  =  0. 

- BO  \b7  N=1.JH - 

SPr*SRF+CMU<  N )*SJ 
1*?  SJ=SJ+OS 

SRF=SRF/CMF 

*PR*  ( rMA*SR-CMF*SnF  1  /CMR 
nSF*SR-SOF 

- r>5'Br3ffff-5E - 

U( 1 )=V(  1 )+OSF»CTH _ 

U(2)=V(?)-0SF*STH 

U(3)=V(3) _ 

u  <  4 1  «v  ( n-nsR*CTH 

U(Msy<2)+nSR»STH _ 

U( 4 1 *V( 3 ) 

00  151  1  =  1.3 

- CJTfeTTWTTTn - 

U(  9+ I )  =V( 3  +  1  ) 

U(  12+1 ) =V ( 6+ 1 ) 

151  m  14+1  1  =V  t  4  +  I  ) _ _ 

II  (  1  Q 1 «V ( 1 )+(5R-SH)*fTH-4F( JHl+STH 
U(?0 )«V(?l-(SR-5Hl*STH-nF{ JHI*rTH 

IT?  ?i  1=1/(101 

RF=(SRF-SH)*(SPF-SH)+OF< JH)*DF( JH) 
PR=(S8R-SH)»<SBR-SH1+DF( JH)*OF( JH) 
SJ=DS/2. 

- nrrc(- - - 

DO  15?  N« 1 » JH 
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C  IF  =  CI F+CIM1N )+CMU(N)*SJ*SJ 

153  SJ=SJ+DS _ 

CIF=CIF-CMF#SRF*SBF 

ClR  =  Cl/>-ClF-CMF»(  SR-SRF  l*(  SP-SPF  )-CMR»  (  SR-SPR  )»(5R-SRR  ) 
G1=(ClF/CMF+rLM*RF) /TM 

Q6=(CIR/CMR+TM»RR  )  /TM _ 

nXPrUm-U(lR) 

OYF=U(  ?0  HH  7  > _ _ 

MRTN=2  .  . . 

nTH=0, _ 

16n  pr>TM  =  DTH 

PFMrFN _ 

r>0  161  1  =  1,21 

161  PUjUfUj  1 1 _ 

no  162  1=1,6 

U(6+I 1=U(6+I  )  +  (U(  12±j  1+PU(  12±J  I  I  *nT/  ?  . _ 

162  Util  =  PU  (  I  )  +  <U(6  +  I)  +  PU(6+m*DT/?, 

THT =U ( 3  ) _ 

no  163  1=1,9 

1  6-»  H(n=A(  I  ,1  )  +THT»  (  A  (  1,21  +  THT »  (  A  (  T  ,  3  )  »T>  T*A  U  ,4  )  )  I _ 

STH=SINF(ThT) 

rTM=CQSF(TMT  1 _ 

IP  (IHT)  167,168. 16  8 

167  M ( 1 ) =SH _ 

158  T1=SB-1./H(1 ) 

H(41=l./H(4) _ 

T2=H(2)*STH-T1*CTH 

m(21=H(?1*CTH+T1»STH _ 

H( 1 )=T2 

5THR  =  S1NF(U(6  )  ) _ 

rTwP  =  rOSF(U(6  )  ) 

U(?1 )*H(?»-nSF»STH-U<2) _ 

Z=U(21  ) 

ZO= ( U ( 2 1 ) -PU ( 2 1 ) 1/QT _ 

DO  164  1=1,2 

U(  Io)=PU(19)+(U(7)+PU(7)+(U(o)+PU(9)  )*DYF)»DT^2. _ _ 

U(90  1=PU(?0)  +  (U(81+PU(R)  +  (U(  0  1+PU1P) l*nXF)#DT/2, 

nXP=U( 1 1 -U (IQ) _ 

RYP=U(70)-U(?1 

DXR=U( 191-U14  ) _ 

164  DYR=U(5)-U<20) 

V<  1 )  f  uni _ 

V ( 2 1  *  U(2  I 

V 1 3 1 SU ( 3 ) _ 

V(4)=U(7) 

V(*l  =  U  <  8  ) _ 

V  (  6 1  «  U  (  9  ) 

T=T+DT 

- £6  TO  ffl - 

165  DTH*U.(  6  /-U  (  3  ) _ 
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JR(  ABSR(nTH)-rr>TH1  IMS*  1  f,f*  1  A7 

147  PRINT  14A _ _ 

1AA  RORMAT  ( 1  AH  FUSFLAGF  RR F A<  ) 

_ GO  TO  70C _ _ _ 

144  DMMrSTRR*ZW(  JH)«(  .4-11.1  *OTH*OTH  ) 

TpnTH-PnTm  1  AO  >17  1,171 _ 

14«  Phm  =  _rpm*PhVI 

171  G4=(0XR*0XR+0YF*0YR) _ _ _ _ 

0?sO£»*RM 

F 1 =U( 9 ) *U( 9 ) _ 

F?  =  U<'l?)*U<12) 

F-a-nxR*  (  H*OYF  +  F?*nYR  )  -DYR*  (  F  1  »PXF+F2*0XR  ) _ 

GA=PHM/ (GMF*TM) +FN# (U(  10 1 / ( rMR*TM)+r»XR/rMR  )+FT*  (  (U( 20  )+. A*Z , / 
1  ( GMC*  TM ) -nyF /GMP 1 +RM*R  A 
ntys-pH*A/  (rMR#TM)-FM#nXF/rMF  +  FT*nyR/GMF+F7 

r>FN  =  G1  »Q4-G4»G7  _ _  ___ _ 

U(lMs(GA*r,5-G?*0M/r>FN 
U ( lA)=(G1*G6-GT*G4)/nFN 

F4rRl#PYF+F2*DYR 

R4=nYF*U(1A) 

F7  =  nYR«U(l A  ) _ 

FA=nXR*U ( 1  A  I 
Ff>-nXR*U<  1  A  ) 

U(  I  7  i  =-  (FT+CMR*F4  )  /CMA-ri.M#  (  Fft  +  F7  ) 

UI161 =-( FT-CMF*F4 ) /CMA+TM# I  F6+F7  ) 

U ( 1 4 )  =  ( FN+GMR*FA ) /CMA-CLM*  <  F8  +  F9  ) 

U< 17)=(FN-CMF*FA)/rMA+TM*{F8+F9) 

5J=PS/2. 

no  17A  J* 1 *NS 

IF  (J-JH1  174.174. 1?A 

1^4  tr  (  j  i  =-U(  1  a  )*CTM+U(  14)»5TH+U<q)*U(Qi«(SRR-SJ) _ 

—  tQ  J71 

174  A0< J)=-U(16)#CTHR+U( 17)*STHR+U( 12)*U( 12)«(SRR-SJ) 

17A  SJ=SJ+DS 
GO  TO  1  40 

700  PRINT  ?  0 A  *  NRUN 

20A  FORMAT ( 1H1 »20X.  1  4H  INPUT  DATA  RUN. I*//) 

PR  t  NT  704.VEl.SNK.THTP.FMU 

706  FORMA  T  (  AH  VRL.F7.7.AX.4H  SNK.F4.7. _ 

1AX.4H  THT.F6.4.4X* 11H  FRirT  CORF, FA.?/) 

PRINT  20S.STRB.CDTH.BETA _ 

704  f6rMAT(11H  CR  STRFSS  .F10.A.6X.7H  CR  ANG.F6.A.6X. 

1  l  2H  DAMP  FACTOR.F7.3/  ) _ 

PRINT  206.H3P .H4P.H5P 

204  F0RMATI8H  FN1  «ZS/ . E  10.4  .  I  OX  .  1 1 H  FN2=ZG»ZG»  .F 1  0.4  . 1  OX . _ 

116H  FN3»ZG*ZGD*ZGP/.F10.4/) 

PRINT  707.H6P.H7P.H8P  _ 

?07  FORMAT  ( 7H  FNG*1 /.F10.6. 10X.11H  FT  1 =XD#ZG* . FI  0.4 . 1  OX . 

114H  FT2*XD»XD»ZG».E10.4// ) _ 
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PRINT  208 

?08  FORMAT  ( 20X  » 1 5H  DEVcL0PFP  DATA//) _ 

PRINT  ?09. ( V<  I  )  .  1  =  1 .6) 

209  FORMAT!^  X  *F7.2»5X*3H  Y  *Ffe*?»^X*SH  T HFT A t F6 . 2 » 5 X 5H  XDQT.F7.2* 
15X.SH  YDOT.F7.2.5X.9H  THETADOT . F6 . 2 / ) 

ZS=ZSP _ 

GL  =  V  (  D-H1P 

PRINT  410*ZG*ZS»ZGM*GL*STRM _ 

410  FORMAT! 12H  GRND  DEFORM*  F5 .2 ♦ 5X  »9H  FUS  DEFl , F5 . 2 . 5X * 

11  MAX  GRND  DFcRM<F6»7*5X»7H  or  poyr  ,  F  7 ,  ;  ,  c,  x  *  1  OH  MAXSTRFSS  *F9. 2/  1 
PRINT  409»  (QMX(  I)  1 1  =  1  tl»  )  .  JST 

409  FORMAT ( 29H  GENERALIZED  COORDINATES  Q  1  =  .  F 1 0. 5 . 3X * 4H  O?  =  , F I  0 . 5 . 3X * 
14H  Q3=*F10.5»3X.4H  04=  »  FI  0.5  »3X  •  1  2H  MAX  STP  100*1*?) 

STRL  =  ARSF(  STRL) _ _ 

PRINT  192»STRL*JSTL 

192  FORMAT  (10X*16H  MAX  TENS  STRESS*  F9 . 2  *  3X  >  1  6H  LOCATION  NUMRFR  *15  1 
PRINT  40 

40  FORMAT  (lHl.llVl  INPUT  DATA) _ 

PRINT  145* ( ( A (  I  .  J)  .  J=1.4) *  1  =  1 *9 ) 

PRINT  145*  FC « SC _ 

PRINT  145*  CMA*CI  A.SP.STRR.CDTH.RPM.CNS.DS.FMU.DT 

NRUN=NRUN+1 _ _ _ _ 

REAP  10.VEL.SNK.THT 

I F ( EOF  *  60 )  999*  11 _ 

999  STOP 

FND _ 


INPUT  DATA  FOR  CRASH  PULSE  SIMULATOR 


Aircraft  -  CV-7  (Operational  Light) 

VEL  ■  120.0 

SNK  =  20.  0 

THT  =  0.  2094 

THTD  =  0. 

A(I,  J): 

A(I,  J) 


1 

2 

3 

4 

1 

. 03704 

3. 1 196 

-8. 557 

11. 005 

2 

4.  64 

-  1. 672 

0.  0 

0.0 

3 

8.  0  x  10“7 

0.  0 

8.  0  x 

10 

-6 

0.  0 

4 

1. 924  x  10"7 

0.  0 

33.  6  x 

10 

-6 

0.  0 

5 

-4 

8.  0  x  10 

0.  0 

135.  0  x 

10 

-4 

0.  0 

6 

2.  0  x  10'5 

0.  0 

42.  2  x 

10 

-5 

0.  0 

7 

0.  0 

0.  0 

0.  0 

0.  0 

8 

40000. 

2400. 

0.  0 

0.  0 

9 

20.  0 

0.  0 

0.  0 

0.  0 

CMA  ■  900.0 

CDTH 

= 

0.  30 

CIA 

5 

=  1.  53  x  10 

RPM 

5 

Y° 

SB 

a  27.  0 

CNS 

= 

26.  0 

7 

STRB  =  7.  0  x  10 

BETA 

40.  0 
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DS  =  26. 

FMU  =  0.  75 

DT  =  .001 

OUTPUT  OF  PROGRAM  VIBRAT  (REFER  TO  APPENDIX  VI) 
CMU,  CIM,  AS,  DF: 


CMU 

CIM 

ZM 

DF 

1 

9.0 

200 

.050 

2.  0 

2 

24.  5 

250 

.  050 

3.0 

3 

28.  5 

480 

.050 

4.  0 

4 

24.5 

390 

.0608 

4.  5 

5 

14.0 

250 

.  1544 

4.  5 

6 

8.9 

250 

.  1831 

4.  5 

7 

19.0 

254 

.  1821 

4.  5 

8 

211.  3 

340 

.  1731 

4.  5 

9 

210.3 

250 

00 

00 

4.  5 

10 

208.  0 

198 

.  252 

4.  5 

1 1 

18.  5 

240 

.  2375 

4.  5 

12 

19.  5 

380 

.  2282 

4.  5 

13 

22.  0 

413 

.  242 

4.  5 

14 

10.0 

320 

.  1748 

4.  5 

15 

8.  0 

310 

.  180 

3.  3 

16 

4.  0 

110 

.  185 

3.  1 

17 

4.0 

65 

.  150 

3.  0 

18 

1.5 

48 

.  1168 

3.  0 

19 

3.0 

48 

.112 

2.  5 

20 

4.5 

27 

.  092 

2.  5 

21 

4.  5 

20 

.  0767 

2.  5 

22 

19.0 

20 

.  0767 

2.  5 

23 

19.5 

950 

.60 

2.  5 

24 

4.  5 

3500 

.  520 

2.  5 

25 

1.0 

3500 

.  10 

2.  5 

26 

1.  25 

3 

.  10 

2.  5 
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Developed  Data  For  Crash  Pulse  Simulator 


Final  coordinates:  (at  rebound) 
X  =  -2.  37  ft 
Y  =  7.  45  ft 

9  =  0.  13  rad 


Final  velocities  and  accelerations: 
X  =  78.  59  ft/sec 
Y  =  -2.  17  ft/sec 
0  =  -0. 50  rad/sec 


Fuselage  deflection 
Maximum  ground  deformation 
Length  of  groove  formed  between 
contact  and  rebound 


0.  78  ft 
0.  45  ft 

20.  8  ft 
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APPENDIX  V 

CARGO  SIMULATOR  FORTRAN  LISTING 


PROGRAM  CRASH 

n  tension  vnoitPvi  io)*a(p*4)  *h(9)  .worn  30  >  .01201  * _ 

1PO(>0)  *PH7  (32  )  «PH7(37)  .PH*  (  V  )  *PH4 (32)  *CMU(3?  I  *ECA( 32  >  *FCB(  3?  ) 
inr(3?)*rTM(3?)*Ar(37l*U(?n*PU(7l)  »OM(  4  )  *FM(  a  )  *SNn(  A  >  »CSD( 4  )  » 

1  rMOMU  ) .PLOT ( 71  l*°T(4|.0MXU).ZM(V) 

10  FORMAT ( 6F10.6  ) _ 

RFAn  10*  v-;l  »snk*  tht  *thtd 

_ RFfln  1C*  ( ( A ( I  * J  )  » J  =  l*4 ) »  1  =  1  *9  ) _ 

RFAn  10*  CMA.CI  A  .  i>R  .  STRB  *COTH  .  RPM  *CN  S  *  DS  .  FMU  *DT 

PFflr>  )0,CMU«CIM,DF>PHl*PH7*PH3*PH4.FCA*EfB*ZM _ 

PFAH  10.OM.CM.RFTA  “ 

PFAn  300* (PLOT (  J) * J=1 *71 1 * RLNK *PR^*flST  *  FO*  PLS _ 

300  FORMAT  (71A1.5A11 

R FAQ  1  0  *CLN  *  RLL  *  FLL _ _ 

RF AD  9.MNTH.NDAY .NYR.NRUN 
9  FORMAT  ( A  I  5 ) 

- 5MIM.5 - 

S r=.7B _ 

11  no  i 2  * J=1 *  1 0 

12  V(J>=0. _ 

PU 

FLL=FLL*32«  2 _ 

LN=rLN 

V  (  3 ) »THT _ 

V ( 4 ) =SGRTF < VEL*VFL-SNK*SNK ) 

Vm=-SNK _ 

V(6)«ThTD 

_ no  14  i«i.? _ 

14  v<n  =  A(I.l  )+T'1T*(A(1.?1+ThT*(A(  !.»3»+THT*AU*4m 

FN  =  9. _ 

IF  (THT)  70*71.71 

-40-V44 WW44t^U4 *  1  >-600.«ThT  ) - * 

71  Tl=SB-l./V( 1) 

T7=V(7)»SINF(THT)-T1»C0SF(THT) _ 

V ( 7 1 =V ( 2 1  *f OSF (THT)+T1*SINF(THT) 

V(  11=T7 _ 

FT=FN  . ’ 

MXT=  1 _ 

NCR=  1 
MRTNs 1 

- VhtP>  »  ThT - 

no  17  IQ.9 _ 

17  HIP  =  A( T *1 )+THT*( A ( ! .7 )+THT* ( A (| .3 )+THT*A ( I .4) ) 1 

H1P  =  H(3)/(1.  +  4.»THT»THT) _ 

H4P=  ( 1 ,+4.*THT*THT ) /H( 4  I 
H5P  *  H  (  5 1 

- MAP=Hl  Ai  /  ( 1  .+'8.^THT*THT » - 

H7P*H  (  8  1 
H8P  =  H( 9  ) 

H1P=V( 1 ) 
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&nmfpm  v  aWiffi-w  mi  wsafei 


DO  13  J  »  l.?0 

13  Q<  J ) *0. _ 

T»0. 

NS«CNS _ 

0«MI  *  0M(3|  -  BFTA  *  PPTA/4.0 

OM ( 4 )  8  OM ( 4 1  ~  RPTA  »  3FTA/4.0 _ 

HO  15  I  «  1.4 

OMX ( 1 )»Q. _ 

CMOM ( I ) »SQRTF<OM< I  1 )*CM( I ) 

SNn( 1 ) a S I NF ( SORT F  ( OM  ( I ) }»QT  ) _ 

15  CSr>(  !)«COSF(SORTF(OM1  m*r*T) 

5WCHg  -1.0 _ 

no  -»oi  j  =>  i.ns 

*r< ji  «  o. _ 

no  16  !  *  3.4 

SNO  < I )  -  FXPF( -BFTA/2.  »  DT )  »  SNDll) _ 

16  CSn(l)  *  PXPF1-RFTA/2.  *  nT)  *  CSWI) 

_ ?T<H«0, _ _ _ 

RT(?)«0. 

PTC*)«BFTA _ 

RT  ( 4  ) *RFTA 

ZG  =  0. _ 

FNC  *  0.0 

FCg5.0 _ 

FU  =  2.0 

SH»40. _ 

ZDG  =  0. 

7GM  »  0. _ 

STRM  *  0. 

PRINT  700.HNTH.NOAY.NYR.NRUN _ 

200  PORMATC 1H1  .20X.6H  OATF  . 2 14 ♦ ?0X . 6H  RUN  .15//) 
PRINT  201 

201  FORMAT  (5X.31H  TIME  CARGO  ACEL  FLOOR  ACEL.3X. 

12RHC)  FLOOR.  Ul  CARGO) _ 

MSW  *  1 

PVR  ■  0. _ 

VR  *  0. 

PAR  »  C. _ 

SR  *  0. 

SRM  =  0. _ 

SRC  *  0. 

STRL«0. _ 

AR*0. 

MCR»1 _ 

00  1 12  1*1.3 

1 12  v(  n»v(3+i  i*r>T+v(  1 1 _ 

1 10  PFNsFN  U 

DO  111  I»1.10 _ 

111  PV( I )«V( I) 

THT  «V( 3 1 
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no  m  1  =  1*9 

115  Htn=A(I.1)+THT»(A(!,2)  +ThT»(  A  <  1  *  3  )  4-ThT  »A  t I  .4 ) ) ) 
IF  (ThT)  117*118.118 

117  H ( 1)=SH 

US  TiVSffi\T7HTTi - 

MU)  =  1,/H(4  ) _ 

T2=H<  2  )*SINF ( THT  I-T1*C05F( THT 1 

H(7>=H(7)»rOSF(THT)+n«MNF(THT) _ 

w( 1  )  =  T? 

V( 10 ) gH ( 2) -V ( 2  ) _ 

ZO=(VI  10  1-pS/  i  id  )  )/DT 

Z  =V ( 10) _ 

127  IF(T-.8)  180*700.700 

18C  I F ( V ( 4  )  )  114*114.181 _ 

181  I F ( Z  )  114.114.182 

18?  IF(V(5)-.R)  1 18.700. 700 _ 

1 14  PRINT  lift 

118  FORMAT (18H  REBOUND  COMPLETE  ) _ _ 

GO  TO  700 

118  SWCH  =  -1.0  »  SWCH _ 

IF(flCLN).  870.571  .871 

570  ACLN=0. _ 

571  IF(VR)  500,501.802 

501  TF  (AC(LN))  505,^04.504 _ 

804  IF  (SRM  -  SRI  505,505,508 

505  IF  (AC(LN)  -  RLL )  600.600,507 

507  AR  s  ACILNI  -  RLL 

508  IF  (VR)  517.518, 519 

5TT  aR  =  AR  +  SMU  »  aclN - 

GO  TO  520 

519  AR  =  AR  -  SMU  *  AClN 
GO  TO  520 

518  IF  (ArSfIaRI  -  SMU  *  ACLN)  521,521.522 

521  AR  =  0 

- T0  yT-4 - 

522  AR  =  AR  *  (1.  -  SMU  *  ACLN/ ABSF ( AR JJ 

550  VR  =  PVR  +  I  AR  ♦  PaR)  *  DT/2. - - - 

1 F  (VR  *  PVR)  509.510.510 

509  vfc  =  0. - 

AR  =  0. 

- <*,c  Term - 

510  SR  =  SR  +  (VR  +  PVR)  *  DT/2. 

IF  (SR  -  SRM)  511.511,512 

512  SRM  =  SR 

511  IF  (SR  --  R L  f  ~~~ V  I  •  5  1  4 
518  SRL  =  SR 

5 14  PVR  =  VR 
PAR  *  Ai 
GO  TO  600' 

506  AR  =  AC(LN) 


GO  TO  608 

603  I  F  (SR  -  SRL  )  516.615.506 _ 

616  Tp  (_  ai-ILW)  -  FLU  #>00 . 600 * 6  1  6 

516  AR  =  AC(LN)  +  FLL _ 

GO  TO  508 

507  IF  ( SRM  -  SR)  507.607.506 _ 

500  JF  (SR  -  SRI.)  516*516.5  06 

600  I F  (SWCM)  51fl.55C.651 _ 

500  FTOP 

571  r,R  =  (  AC  (  LN  )  -  AR  1  Z32..2 _ 

K  =  GO  +  76.6 
IF  <K)  505.505.706 

506  IF  ( K-7 1 )  507.507.505 

507  PLOT  ( K )  =  EO 

5 J 6  GC  *  AC (LN 1/32.2 

L  •  GC  t  76.5 _ 

IF  (L)  5rfl.5CB.50q 

50g  IF  (L-711  510.710, 7Qfl _ _ 

510  PLOT  (L)  =  AST 

308  PLOT  (26)  «  PRD _ 

PRINT  5 1 1 . T • GB • GC • PLOT 

311  FORMAT  (5X.F5.3.2F12.2. 7IflI) _ 

312  IF  IK)  315,315.316 

516  IF  (<-71 )  517.317.715 _ 

70  PLOT  U1  =  RI.NK 

316  IF  (L)  318.318.319 _ 

319  IF  (L-71)  320.370.518 

320  PLOT  (LI  =  BLNK _ 

318  GO  TOi 523.524) »MCR 

523  PZ-ZD«0T _ 

COFsl ,+A.*THT*THT 

H(7l=H(3)/rOF _ 

M(4)=HU)*rOF 

H(6lsH(61/<2.»CQF-1.1 _ 

I F ( ZG  -  ZGM 1  41 1,411,412 

4  12  Z  GM  =  Z  G _ 

ZSP=Z-ZG 

41 1  PZG  =  ZG _ 

IF  (DZ)  119.120.120 

1  20  GO  TO  (800.801  )  .  MSW _ 

800  FUH=H ( 4 ) +20G  *  Z0G/H(5) 

CKs( ( 1.0-SC)»H(81»V(41-H  .0/H(3) )/(?,0  »  FUm) 
ZG=S0RTF(Z/(H(7)*FUH1+CtC*0<1-0K 
ZS  =  Z  -  ZG 
FN  =  ZS/H( 3 ) 

IF  IFN-1./H(6M  803,803.802 

MSU  =  2 - 

801  FUH  «  H ( 4 1  +  ZOG  »  ZQG  /H(5) _ 

IF  (FN  -  1 «/H(6 1 )  804.805, 80 6 

804  FCH  =  H ( 3  1 
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GO  TO  806 

AQ6  pru  a  PC  »  HH| _ 

80  A  rizr.  =  nz/  (  (  ?.#FUH*Z0+<  1  .-SC ) *H <  8  t *V ( A  I  )«FCH+1  .0) 

pqq  zg=zg+ozg _ 

nzs  ='  nz  -  ozg 

DFN  =  DZS/FCH _ 

FN=PFN+OFN 

GO  TO  803 _ 

1  IQ  OZG=DZ*ZG/ (Z-DZ) 

FCH  =  H(3) _ 

GO  TO  899 

flfV*  FNSr=( 1 .-SC)«H(8)»V(A)»2G _ 

IF  (FN-FNSC)  808*121*1?! 

808  FN=cNSC _ 

121  FT  =  FMU  *  FN  +  ( H( 8 )  *  VIA)  +  H(O)  *  VIA)  *  VIA))  »  ZG  *  SC 
ZOG  =  (ZG  -  PZGI/DT 
fTh=CA$F I  V ( 3  )  ) 

sth=s;nf( VII ) ) _ _ 

P=  cN*rTW-FT«STH 
r  =  FT  *  r TH  +  FN  *  STh 
1 2*  GO  TO  I  1 2A.1 AS) *MRTN 
12A  V(7)=-FT/C^A 

- --fleR - 

V  I  8 ) =  FN/CMA 

V(9)  =  (fN*H< 1  )+Fl*(Vi2)+ZG/2. ) ) /C I  A 
00  130  1  =  1  ,3 

VI3+I  )=V(3  +  I  )  +  ( VI 6+  1 1  +PV ( <S  + 1  )  )*oT/?, 

190  V  I T l = V I  I )  +  I  V  I  3  + 1 I+PVI3+I  1  1 *0T /?. 

- - 

00  131  J= 1 *NS 

WOnl J ) =V( 7  )*STH+V( R )«CTH-V  I  A ) *  I Sfl-S J ) 

131  SJ=SJ+DS 

00  132  I  =  1.20 
1  •»?  PO I  I  )  =0  I  I  ) 

- V{ ?  )+2<V?.  )*5Th+hI  i  pkTTh - 

CMM=C*(  I  VI 2  )+ZG/2. ) #rTH+H(  1  )#STH) 

J=SP/0S+1 .0 

Q  < 1 )  =  p  ♦  PHI  I J )  +  CMM  *  (PHI (J  +  l  ) -PHI ( J ) )/PS _ 

0(2)  S  P*  PH2IJ)  +  CMM  *  (PH2I J+1)-PH2(J))/DS 

0(3) aC*PH3 I J  ) _ 

0(A)=C*PHA ( J) 

ACFL  =  V ( 8  )  »  STH  -  V I  7 )  »  CTH _ 

SJ  =  S*  -  OS/2. 

THSO=V(M»VIA) _ 

00  133  J=1 »NS 

0(1)  =  0( 1  )-CMU( J)*WDO(  J)»PH1 1 J) 

0(2)  =  QI2  l-CMUl J)*WDD(  J)*PH2( J) 

Q  ( 3 )  =  0(3)  -  CMU(J)  *  (ACEL  ♦  THSQ  *  SJ)  »  PH3(J) _ 

0(A)  =  0(A)  -  CMU(J)  *  (ACFL  +  THSO  *  SJ)  *  PHA(J) 

133  SJ  =  SJ  -  OS 
00  21?  I  =  1 » A 


j 


5 


lF(T-0T/2.)  210*210*21 1 _ 

210  0(4  +  I)  *  0111  *  DT  *  DT/<6.  *  CM(m 

qua  ♦  n  •  o(4  +  n  »3./( sortfiomi i n  »  dti _ 

GO  TO  213  u. 

211  OT  *  PQ(  16  +  n  »  PQ  (  n  /CMOM(  I)»r>T _ 

0(4  +  n  =  PO ( 4  +  II  *  050(1)  +  OT  *  SND(I) 

0 ( 1 A  ♦  I)  ■  QT  «  CSO(l)  -  PO ( 4  +  I)  »  SND(I) _ 

213  0(12+11=0(1  I  /CM  (  I  )-(0M(  I  )+RT<  I  )*RT(  I  )  /4 .  )  *0 ( 4+  I  ) -RT  (  I  )*PQ(8+I  ) 

I F ( AB5F ( Q ( 4+ 1 ) ) -QMX ( I))  212>212»214  _ 

214  OMX ( I)«AOSF(0(4+I ) ) 

212  0(R  +  I)  =  PQIfl  +  I)  +  (0(12  +  1)  +  PQ ( 1 2  +  1))  »  QT/2.0 _ 

PSTR*0. 

PPSTR.Q. _ _ 

OMNT  =  CMA  -  CMU(l)/2.0 

SJ  s  SR  -  OS/2. _ 

PSHR*0. 

RMM=-CMM _ 

NPP=1 

no  400  J  =  1  ,N5 _ _ 

AC( J)=ACFL  +  THSQ*SJ+0( 15  )*PH3( J)+0( 16 )*PP4( J ) 

0 n  TO  (143.144) .  MXT _ 

141  A0N*V(7 )*STH+V( R 1*0TH-V(0>#5J 

IQ  ( a  -  LN)  146» 147  « 146 _ 

147  ACLN  =  ACN+32.2*CTH 

146  SHR=PSHR+CMU( J)»ACN _ 

GO  T0(40l .402) *NPP 

401  IF  ( SP-SR+S J-DS/2 . )  403.403.402 _ 

403  NPP*2 

SMR.SHR-P _ _ 

40^  BMM=RMM-(PSHR+SWR)»0S/2.+riM( J)«V(P) 

PSHR=SHR _ 

STR=ECA( JJ «Q(5 l+ECR ( J)*0(6) 

STRN=<3.»STR+6.*PSTR-PPSTR)/8. 

IPfSTRL-STRN)  101.191.190 

IRQ  STRL=STRN _ 

JSTl *J 

1R1  I F ( STRN-ST  RM )  413.414.414 _ 

414  STRM  «  STRN 
JST=J 

413  I F ( STRN  -  STRB )  140*141.141 

140  PPSTR»PSTR _ 

PSTR=STR 

GO  TO  144 _ 

141  JH=J-1 

SH  =  SR  -  SJ  -OS/2.0 _ 

PRINT  142.JH 

142  FORMAT ( 2 1H  YIELD  HINGE  AT  STA.  .13) 

- MTTO - 

144  CMNT  *  CMNT  -'CMU(J)  +  CMU(J  +  l)»/2.0 

4H0  5J  =  5J  -  DS -  -  - - 

145  FORMAT  ( 1X.10E11.4/I 


- T=T+I!)T - 

GO  TO  (  1 10  1 1 49 ) »MXT 

T7~  d’MFrir. - 

T  =  T-DT _ 

00  150  N=1|JH 
] 50  rWF=CMF+rMU(N ) 
rwRrCMA-CMF 
CLM=CMR/CMA 

- TM=cHr>(‘MA - 

RM=fMR/CMF 

Sj=oS/2. 

SPP*0. 


on  15  2  N  =  1 ♦  JH _ 

c,bc-  =  Srf+CMU(N)*SJ 

152  S J=SJ+DS _ 

SRF*S0F/CMF 

SRR=(CMA»SB-CMF»SPF)/CMR _ 

O.SF  =  SB-SRF 

ORR=SRR-SR _ 

IK  1  1«V( 1 )+0SF*CTH 

U  (  ? )  =  V  ( 2  )-OSF*STH _ 

um«v<?)  "  “  ~  . . 

U(4)=V( 1 )-DSR»CTH _ 

U(5)=V(2)+DSR*STH 
U  (  6  >  XV ( 3 ) 

- b*  1 51-^173 - 

U ( 6* I  1 gV( 3+1  ) _ 

u(o+n*vi3+i  i 

U(  l?+l  )=V(4+I  ) _ 

1R1  U(  l«  +  n«V(A  +  !  | 

U(  1Q)«V(  11  +  (S8-SH)#CTH-r>F<  JH)»$TH 

- Ul  J0)«VI  J)-(5B-5H)*5TH-nF  (  jh)*cTh - 

U( 21 ) = V ( 10) 

RF=(5RF-SM)*(SRF-SH)+OF( JH1+OFI JH I 
RR= ( SBR-SH ) * ( SRR-SH 1 +DF  «  JH  »  *OF ( JH) 

SJ=o5/2. 
n I f=0. 

00  153  N=1 i JH 

rTP=ClF-»C1M(N)»CMU(N)«SJ»5J _ 

153  SJ=SJ+OS 

C  T  P=C I P-CMF»SRF*SBF 

C IR-C I A-CFF-CMF*!  SR-SBF )  *  (  SR-SBF 1-CMR*  ( SR-5RR 1* ( SB-SBR  1 
G 1 = < C I F/CMF+CLM«RF ) / TM 
GR= (CIR/CMR+TM*RR) /TM  • 

O  X  F  «U ( 1 )  -U ( 1  9 ) 

OYF=U ( 20 )-U ( 2  1 
MRTN=2 

- F5TH57T; - 

160  POTH  «  OTH 

- PFNsrw - ; - 


■\ 


\ 
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161  pum*u(it 

_ PO  162  1*1.6 _ 

u(6+n  *pu<  6+i  >+pu<  i2+n  *dt 


16*  H(  n*A(  I  *1  )  +  THT»(  M  I  »?1+TmT»(  A  (  I  .1 )  +  fHT*A  (  !  •4))) 

STH*SINF(THT) _ 

CTH*COSF(THT) 

IF  (IHI 1  157.158.168 _ 

1 67  M ( 1 ) =SH 

168  T 1  =  SB-1 . /H (  1  ) _ 

H(41*l./H(4) 

T7sH(?1»STh-T1*CTH _ 

H(2)*H(2)»fTH+Tl*STH 

H ( 1 ) «T2 _ 

STHR  =  SINF (U46  )  ) 

CThR=C05F(U(6 1  ) _ 

U(21 )«H(2)-nSF*STH-UI2) 

Z*U(21 1 _ _ _ _ _ _ 


DYF=U( 20 )-U( 2 1 
DXR=U ( 1 9 ) -U ( 4 ) 
164  HYR=U  ( 6 ) -U ( 20  1 
V( 11  =  U(1  ) 


V(2)  *  U<2) 


167  PRINT  168 

168  FORMAT  (15H  FUSELAGE  BREAK) 
GO  TO  700 


166  PHH«STRP''ZM(  JH  )  »  (  »5~1 1  »  1*HTH*PTH  ) _ 

IF(OTH-PDTH)  16^* 171*171 

169  PHM=-RPM»PHM _ 

171  G4*(DXF*DXR+DYF*0YR) 

G2«G4*RM 
F 1 =U  (  9 ) *U ( 9  I 
F2=U(  12 )*U( 12  ) 

— Fs=^xH»iFr*DVF-t-rg«c>vK)-&vR#iri»6yr>p^'»c>XR) - — 

G3*PHM/(CMF»TM)+FN»(U(  19 )  /  ( CMF»TM  l+DXF/CMF  )  +FT»  (  ( U  (  20  )  +.  5»Z  )  / 


1 <CMF*TM)-DYF/CMF )+RM*F3 

G6g-PHM/(CMR«TM)-FN«DXF/CMF+FT»0YR/CMF+F3 _ 

0FN=G1*G5-G4*G2 

U!  15,)«!G3»G5-62»G6)  /DEN _ 

U!  18)  =  <G1*G6-G3*G4) /OEN 

F6  =  Fl»r>XF+F2»DXR _ 

F5=F1*DYF+F2#DYR 

_ FftsDYF»U  (15) _ _ 

F7=nYR#U (18) 

FRsnXF»U( 1 5  ) _ 

F9=nXR#U< 18 ) 

U<  13  ’,=- (  FT+CMR»F4  >  /CMA-CLM* ! F6+F7  ) _ 

U(  16  =-(FT-CMF*F4)/CMA  +  TM*(F6  +  F7) 

U(  14)*(FN+CMR»F5)  /CHA-CLM»(  F8  +  F9  ) _ 

U(  1})s'<fN~cMf^5)  /Civ'A+TM*(F8  +  F9) 

SJ-nS/?. _ _ _ 

no  173  J=1.NS 

IF  (J-JH)  174*174*175 _ 

1  74  AT ( J ) =-U (  13  )*CTh+U(  14)#STh+U(9)*U<9)*(SRF-SJ) 

GO  TO  173 _ 

1  7S  AC  (  J)  =-U(  if,  )*CTHR+U(  17)  *SThR+u(  1  2 )*U( 12 ) * (  SPR-SJ  ) 

173  SJsSJ+DS _ 

IF(CLN*3.-SH)  572.572.573 

5  7?  ACLN=U( 13)»STH+!U(  14  )+32 . 2  ) »CTH-U ( 1 5  )  » (  SFF-'!LN»3 .  ) _ 

GO  TO  160 

5  73  Ar|_N=U<  1 6 )  *STH+  !  U  (  1  7  )  +3  ?  .  ?  )  #CTH-U  (  18  )  »  <  S(3R-CLN»3  .  ) 

- WTO  UO - - 

700  AC(LN)=0. 

T=T+DT 
MCR*2 
AOLN*32.2 

GO  TO  115 _ 

524  T  =T+DT 
IF(VK)  525*525*  115 

525  PRINT  203 » NRUN 

203  FORMAT! 1H1 .20X.15H  INPUT  DATA  RUN. 15//) 

PRINT  204.VEL.SNK.ThYP.FMU 

204  F0RMATI4H  VEL.F7.2 .6X.4H  SNK.F6.2. 

16X.4H  THT.F6.3.6X.11H  FRICT  COEF.F5.2/) 

PRINT  205.STRfl.CDTH.RFTA _ 

205  FORMAT ( 1 1H  CR  STRFSS  .F10.3.6X.7H  CR  ANG.F6.3.6X. 

1 1 ?H  DAMP  FACTOR. F7. 3/) 

PRINT  206.H3P.H4P.H5P 

?06  FORMAT  (  8H  FN1-ZS/.E10.4.10X.11H  FN2=ZG*ZG* .E10.4.10X. 

1 T  6H  PN3  =  ZG*Z"GD*ZGD/ *E10,4/ I 

PRINT  207.H6P.H7P.H6P _ 

207  FORMAT  ( 7H  FNC=  1  / • E 10.4 . 10X . 1 1H  FT1=XD*ZG*.E10.4»10X. 

114H  FT2«XD»XD*ZG*.E10.4//) 

PRINT  200 

?0B  FORMAT! 20X.15H  OFVELOPFD  DATA//) 

- PRlNI  m.JVUi. I»1.6) - 


209  FORMAT ! 3H  X  .F7.2.5X.3H  V  .F6.2.5X.6H  THETA .F6.2 .5X5H  XD0T.F7.2. 


1«X.*H  Yr>0T.F7.2.5X»9H  THETADOT .Ff . 2/ ) 

ZS«ZSP _ 

GL«V<1)-H1P 

PRINT  410»ZG«Z5>ZGM >GL»5TRM _ 

410  FORMAT ( 1  2H  GRNO  DEFORM* F5. 2  * 5X  *9H  FUS  DEFL ♦ F5 * 2 1 5X » 

115H  MAX  GRND  PEFRM t Ffe. ? . 5X » 7H  GROOVE »F7. 2  * 5X » 10H  MAXSTRESS » E9 . ?/ > 
PRINT  4C9*  lOMXfn  #  I«1»4I«JST 

4QQ  FORMAT ( 29H  GFNFRAL I ^FP  COORD  I  NAT  PS  Q1=»F10»5»3X  >4H  Q?= tF10.5  .?X » 
14H  Q3».F10.5.3X«4H  04  =  «FIO.^*?X*1  ?H  MAX  5TR  LOC*m 

STRL«APSF(STRL) _ 

PRINT  192*STRL  * JSTL 

192  FORMAT  I10X*16H  MAX  TENS  STRESS. E9. 2 »3X. 16H  LOCATION  NUMBER  » I  5 ) 

PrTnT  193 

193  FORMAT! 1H0 . 5X . 1 3HRFAR  LOAO  LIM.5X.13H  FWD  LOAD  LIM.4X. _ 

114H  MAX  FWD  DISPL.4X.14H  MAX  AFT  DISPL.4X.9H  SECT  LOC ) 

PRINT  194.RLL.FLL.SRM.SRL.LN _ 

194  F0RMAT(flXtF8.?»l0X.FB.?*l0X*F8.?#10X*F8.2»oX»l3///) 

PRINT  40 _ 

46  FORMAT  (1H1.11H  INPUT  DATA) 

PRINT  145* ( (  A  ( I  »  J  ) »J*1.4).1«1»9) _ 

pRjNT  145,  pc,sc 

PRINT  145.  CMA. CIA. SP.STRB.CDTH. RPM.CNS.PS. EMU. DT _ 

N«UNrNRUN+l 

RCAP  10»  VFL  »SNK. THT  »CLN.RLL»FLL _ 

I F ( FOF • 60 1  999.11 

999  STOP _ 

fnp 
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INPUT  DATA  FOR  CARGO  SIMULATOR 


Aircraft  -  CV-7  (Operational  Light) 

The  input  data  is  the  same  as  that  for  the  criah  simulator,  (refer  to 
Appendix  IV)  plus  the  following: 

CLN  -  13. 

RLL  *  8. 

FLL  =  4. 

Developed  Data  For  Cargo  Simulator 

Same  as  for  crash  simulator,  plus: 

Maximum  forward  displacement  of  cargo 

relative  to  floor  s  0.  33  ft 
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0.100  8.69  10.93 

0.102  8.70  10.23 

0.104 _ 8.73 _ 9.64 

0.106  t . 75  9.16 


•Will' 
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APPENDIX  VI 

PROGRAM  VIBRAT  FORTRAN  LISTING 


_ PyPP-RAV  VT=»3/*T _ 

n  r:cKr, i  or.  ph  i  3? )  »r 1  ( ?? )  < <sr  ( 32  i  ,^4  ( 32 )  »kljt  ( 7 1 )  »c  -  'j(  3? )  w  ( *?  > . 

1  ° '.'(■»?)  ,  rp  (  ■>?  )  ,  AS(  ■>  7  )  _ 

o  r-.T./vj  (  c  r  1  r;  .  4  ) 

A  o.PS.FPS  „^c,  C|,. _ 

o  codvat  (  ?n  0.4.7  T  l  f  ) 

_ oci>n  p,  H  c  t  c  T  ,  pr  ,  r-.i  j _ 

or»r-  7  ,  (  p|_nT  (J)  »J=1  .71  |  »=>|_MX  »  f  TAP  .OPP 

7  foppat  mm  ,yi  1 _ _ 

ts  1  p^'OOOOC.G 

_ FFA=riS _ 

l"'r  =  A .  74  7  (rN.S#PS.) 

S.j  =  PS/7  . _ _ 

no  qoo  j=i  ,(vic  “ “ “ 

pu  (  j  t  =  rOGF  (  rnr«:j  ) 
mo  c.jr<- j  +  n;  ~ 

or)  v  =  Cm 

CJ  =  OS/ 7.  —  “ 

pi  =  f'X'lll  *  5J/4. _ 

P7  =  (=.  «  O’Ull)  -  P'll(7il/R.  ' 

pp  =  pi  »  p.S/4. _ _ _ 

PA  =  f7  «  PS/4 

_ PO  poo  =  ;,mc  _ 

Pid  =  r(.*|j(  j  _  1  ,  *  Sj 

S.J  =  SJ  +  PS. _ 

PCI  =  (PF1  +  C'-J(J)  *  SJ1  /  7 . 

_ c 1  =  FI  +  PCI _ 

P7  =  P7  +  (FI  -  PF1/7.1  *  PS 

p  p  7  =  -  n  .  r.«D  (  ] '  1 /7. _ _ 

C7=F7+PF7 

o~n  r !,  -  Cj  +  (F?  -  PC7/7.)  v-  os _ 

Pp"i  =  (  FMI I  (  M  S  )  *  (  SJ  +  ^F/4.1  )/7* 

FI  =  FI  +  DF 1 _ 

p  7  =  FI  +  (FT  -  PP 1  /  7  * )  *  PS/?. 

pp 7  =  rr-iU(  f.!S  )/?. _ _ _ _ 

P7  =  F?  +  PF? 

P  4  r  C/.  4-  (F?  -  PF?/?,|  »  PS/7. _ 

'*XT  =  ’  - - 

QP1  V(  1  '  =  (  c  .«rK\i  (  1  )  «pu(  j  i-r-’j  i(  7|i!PH|?|  l  /  Q. _ 

OM(1 )  =  V<  1 I  *  PS/4. 

PP  Pv7  J  =  7. MS _ _ _ 

PV=  (C"'U(  J-l  )epu(j_i  |  +CMU  (  J  )  *HH  (  J  )  1/7. 

V( J)  =  V ( J  -  1)  +  PV _ 

9"?  =  Fr!(j  -  1)  +  (v(J)  -  PV/7.)  *  pp  - 

nr  to  ( QP7 .pQ4) ,»xt _ 

op.  pv=(p.<-cr,i)(MS)«Pu(f'S)-r'"'(Ns-''  1 1  /«.  ' 

VF  r  V(KS)  +  PV _ _ _ _ 

o!-c  =  PM  (MS)  +  (V(N-S)  +  PV/7.)  *  PS/2. 

PFM  =  FI  *  F4  -  F7  *  FA _ 

~ =  (P?  *  RMF  -  F4  *  VFl/HFM  - 
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i*?  -  (ra  »  V c  - 


ZJ  = 

rv<  -  \/|Qh|ii  +  01  *  CJ  *  r?i 


no  '  .  -  J  =  1  .NS 
Du(j)  =  PH{J|  +  n  * 


=  :J  +  ns 

-  y 


nPHP  =  <°r(J  -  it/ciu  -  1)  +  R*1  ( J )  /  S !  ( J  )  )  *  n-^/2 
_ pup  =  pup  +  rip  hp 


0  ->  PM(J)  =  HH(J  -  1  1  4-  !PHP  -  0PHP/2.)  *  OS 
_ pprf.'T  oo.pu 


C-0  TO  avi 

on/. 


Ic  ( DO,"  -  PPM  QOft.O08.Rll 

q  is  no  o.iQ  J  =  l  . 


PLOT  ( ■>«.  1  =  °P0 

l-i.T)  sOi.eoj.so. 


put  :1t  ,+PM  (  .11  *?  .  . 
K'  -  PUT  4-  1  .  R 


.PLOT 


1C  FOP'  AT  (  n.4X.c8,4»lJX.7141//) 
Q o  PLOT ( K 1  =  8LN< 


00  TO  ( °?1 .9?? ) »MSW 
Qr>l  VSv.’  =  2 


POM=L . 


ovn  =  rv i  +  r"J(J)  *  ppfji  *  PMfji 

014  PHMJI  ■-  Ph(J1 _ 

PRINT  14.  CM1 

PRINT  !■» _ 

n  roR»'4T  I  1H1  . 1  ?H  £PfONn  mcof  ) 

_ pJ=nS/?. _ 

rop=o.4?  /  (rNS*nsi 

no  o 1 R  J=1 .ns _ 
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r  1  *•  J  =  1  .MS 

-  rn  4.  fVtllJl  »  pu|  II  *  DHM.lt 


'  =  rw/ri/t 

O  o  t  7  J  =  l  . 


OhI  Jl  =  PH(J)  -  A  *  PHA(J) 

u  (  om  .ot  a  t . m  - 


c«rT=l ./PW( 1 ) 

=  1 . 


Du( J)=pu(J ) #FACT 
cnoyjr  (cx.inn  t.M 


no  TO  oOl 

rM-ii. 


$=.0001 

1  0=0  J  =  i .  - 


ru  =  rM  +  CviU(  J)*PH(  J1*PH(  J1 
PplMT  1 L .  cv 


I.HMs  1 

cnr.rnr/5. 


e J=oS/7 . 
POM=0. 


on  o$o  j=i  ,m.$ 
Pu(.M=ro$r(rOF»$J  ) 


«.J=0J  +  ns 
POr^'T  oo,Ph 


rM!.JM=  (  =  {  1  t  *PH  (  1  1  -<■«!.!(  7  1  (71  I  /P  . 

no  009  J  = 


nrMUM=(rwU(  J-1  )*PM(  J-1  1  +  r*«U  (  J  1  *PH(  J  t  t /9  . 
QNijv.=rr;u,<  +  nr  n  u  m 


nr ’IU(;=  (  =  .*CMU(NS)*PH(\'S  1  -C^UfNS-l  )*PH(NS-l))/8. 
ri = ( cnum  +  dcnumi/f 


DM=1 .  /  (PH( i 1-ri 1 
r-0  Tu  (Q77.Q7H).  M$W 


no  079  j=i ,ms 

PM  (  ,|  1  =  (  PM(  ,n -r  1  t*0'< 


nov=OM-POM 

pptmT  go.rwLiv.ri  .oh,4 .cm.omi 


DpM=OM 
PPT  NT  QO .PH 


PH(1)=  CF(1)*0S  /  U.*P*A?<11) 
no  070  1=7. NO 


nPH=  (CPIJ-11  /  A  S ( J— 1 1  +  CF(J)  /  A  S ( J 1 1  *ns  /  (?.*F) 
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INPUT  FOR  PROGRAM  VI  BRAT 


DS 

3.  0 

EPS 

0.  10 

NS 

26. 

MSW  = 

1 

CMU  = 

refer  to  Appendix  IV 

AS 

El 

EC 

I 

.05 

.  7E8 

4.  33E9 

2 

.  037 

1. 455E8 

5.  76E9 

3 

.  037 

3.  0E8 

5.76E9 

4 

.033 

4.  89E8 

7.  74E9 

5 

.  1190 

15. 20E8 

9.48E8 

6 

.  1122 

17.  5E8 

9.  17E9 

7 

.  1165 

17.  6E8 

9.  28E9 

8 

.1118 

15.8E8 

8.  80E9 

9 

.  1360 

29.  3E8 

9.  97E9 

10 

.  1360 

29.  3E8 

11.  25E9 

11 

.1124 

24.  2E8 

9.  80E9 

12 

.1228 

27.  4E8 

11.  52E9 

13 

.  1110 

23.  0E8 

8.  8E9 

14 

.0971 

15.  8E8 

8.  8E9 

15 

.0784 

8.  85E8 

4.  8E9 

16 

.  0735 

8.  48E8 

4.  4E9 

17 

.  0686 

6.  64E8 

4.  32E9 

18 

.  0660 

5.  15E8 

4.  36E9 

19 

.0610 

4.  13E8 

3.  60E9 

20 

.0680 

3.  21E8 

3.  60E9 

21 

.0631 

2.  82E8 

3.  60E9 

22 

.  20 

2.  82E8 

3.  60E9 

23 

.  20 

22.  1E8 

3.  60E9 

24 

.  20 

19.  2E8 

3.  60E9 

25 

.  20 

!. 77E8 

3.  60E9 

26 

.  20 

1.  49  E8 

3.  60E9 
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OUTPUT  OF  PROGRAM  VIBRAT 


OMU) 

OM(2) 

OM(3) 

OM(4) 

=  3573.  09 

=  21774.42 

=  59345.  91 

=  106076.68 

PHI 

PH  2 

PH3 

PH4 

1 

1.0 

1.0 

1.  o 

1.0 

2 

.8192 

.6066 

.9664 

.9403 

3 

.  6456 

.2713 

.  8666 

.7681 

4 

.  4840 

.0249 

.  7209 

.  5252 

5 

.3359 

-.  1321 

.  6059 

.  3405 

6 

.1991 

-.  2223 

.  5215 

.2123 

7 

.  0737 

-.2535 

.4270 

.0763 

8 

-. 0366 

-.2156 

.  3081 

-.0575 

9 

-. 1299 

-.  1114 

.  1587 

-. 1596 

10 

-. 2082 

.0367 

0063 

-. 2027 

11 

2724 

.2081 

1952 

-. 1961 

12 

-. 3214 

.  3904 

-.3897 

-. 1674 

13 

-. 3550 

.  5713 

-. 5770 

-.  1322 

14 

-. 3715 

.7379 

-. 7762 

-.0876 

15 

-. 3647 

.  8723 

-1.0014 

-.0313 

16 

-. 3273 

.9547 

-1. 2483 

.  0338 

17 

-. 2567 

.9728 

-1. 5021 

.  1031 

18 

-. 1510 

.9177 

-1. 7613 

.  1753 

19 

-.  0079 

.  7820 

-2.  0282 

.  2505 

20 

.  1735 

.5664 

-2.  2802 

.  3221 

21 

.  3909 

.2848 

-2.  5099 

.  3879 

22 

.  6374 

-.0311 

-2.  6511 

.4285 

23 

.9011 

-.3493 

-2.  6904 

.4399 

24 

1. 1703 

-.  6589 

-2. 7065 

.4445 

25 

1. 4424 

-.  9253 

-2.  7116 

.4460 

26 

1.7212 

■  1.0710 

-2.  7136 

.4465 

ECA(J),  ECB(J): 


ECA 

ECB 

1 

5.  165E5 

3.  806E6 

2 

4.  ~?7E6 

2.  760E7 

3 

7.910E6 

4.  016E7 

4 

1.406E7 

6.  318E7 

5 

9.  295E6 

3.  747E7 

6 

1.  132E7 

4.  142E7 

7 

1.  511E7 

5.038E7 

8 

1.979E7 

5.  679E7 

9 

1. 390E7 

2.996E7 

10 

1.715E7 

2.  424E7 

11 

1.805E7 

1.454E7 

12 

1.802E7 

5.  221E6 

13 

1.  618E7 

-2.359E6 

14 

2.  082E7 

-9.837E6 

15 

1.834E7 

-1.292E7 

16 

1.576E7 

-1.361E7 

17 

1.687E7 

-1.620E7 

18 

1.830E7 

-1.830E7 

19 

1.538E7 

-1.504E7 

20 

1.513E7 

-1.266E7 

21 

1.  203E7 

-4.982E6 

22 

7.407E6 

5.394E6 

23 

5.  370E5 

1.  852E6 

24 

3.941E5 

2.989E6 

25 

3.  356E6 

3.741E7 

26 

3.563E6 

4.793E7 
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